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1 . RESULTS OF THE SECOND YEAR RESEARCH 


The report prepared at the end of the first year of this NASA grant indicated work 
had begun on collecting an extensive data base of satellite particle, electric field, image, and 
plasma data. These data were to be used to determine correlations between the fields and 
the particle auroral boundaries. A data base of 8 days of excellent coverage from all 
instruments was completed. Table 1 lists the geomagnetic conditions associated with each 
of the selected data periods. This data set consists of 2 active days, 2 moderately active, 
and 4 geomagnetically rather quiet days. Auroral images were obtained from Dr. John 
Craven at the University of Iowa. A total of 944 images were converted from their raw 
form to auroral emissions in the geomagnetic coordinate system. This data set consisted of 
313 UV images and 63 1 which were a combination of visible images mainly obtained using 
a 577 nm filter. Table 2 indicates how many UV image passes per study day were obtained 
as well as the total number UV images for each day. The intensity, oval size, and fine 
structure track well with the general Kp variations shown in Table 1. Indeed, they also 
track well with the Ae index during each day. As an example, on day 351 the first 15 
hours are extremely quiet and the Ae index enhances to well over 500. During the quiet 
period the images from the first two passes show a faint oval whereas the images from the 
third pass during the enhanced geomagnetic activity are bright, structured, and show an 
expanded oval. 

For each of the days listed in Table 1, both VEFI electric potential data and LAPI 
integrated particle energy fluxes were obtained. On the average between 8 and 1 1 passes of 
useful data per day were obtained. These data were displayed in a format such that either 
the statistical electric field model potential or the statistical precipitation energy flux could be 
superimposed. The Heppner and Maynard [1987] and Hardy et al. [1987] models were 
used for the electric potential and precipitation, respectively. In addition, the auroral image 
intensity along the Dynamics Explorer-2 (DE-2) satellite pass could be computed and 
plotted along with the LAPI precipitation data and Hardy et al. [1987] values. 




An indication of the results of this study are obtained from the Sojka et al. [1991] 
(see attached pre-print) DE data and Time Dependent Ionospheric Model (TDIM) 
comparison. The equivalent LAPI and VEFI comparison with statistical models is given in 
Figures 3 and 4 of this attached pre-print. Figure 3 [Sojka et al., 1991] compares the VEFI 
and Heppner and Maynard electric potential for 5 distinct orbits on day 326. Qualitatively 
the agreement is reasonable. Indeed, this set of VEFI data is in the top 5% of days with 
well behaved (systematic) variations in the set of all DE-2 passes [/. Heppner private 
communication, 1984]. Even so, the quantitative agreement is not good for each pass, i.e., 
for the first three passes the minimum negative potential from the model is at least 30 KV 
less than observed. When comparing the LAPI data with the Hardy et al. [1987] 
precipitation, Figure 4, the differences are even more extreme. The “topology” does not 
show good agreement The model polar cap precipitation is significantly too high, the 
dawn and dusk ovals are too wide, and the maximum intensities are too large. Only the 
latter discrepancy can be adjusted by a rescaling of the absolute fluxes. 

Unfortunately, the other 7 days of data extend this general tendency for the models 
to disagree strongly with the observations and for the variability of one pass to the next to 
be so large that they are different patterns. In summary, even from a visual comparison of 
the 8 days of data the statistical oval and electric potential moiphologies do not look like the 
LAPI and VEFI moiphologies on a pass-by-pass basis. Even adjusting die Hardy oval for 
each pass such that the three hourly Kp was used rather than the daily average from Table 1 
did not produce significant improvement The pass-to-pass differences were still larger 
than tweaking the statistical oval. In the electric field model both the IMF (By and Bz) and 
the Kp indices were adjusted, however, this did not produce a systematic improvement in 
model-data orbit-by-orbit comparisons. For most orbits the disagreement between model 
and data would, at its worst, exceed the value of the model potential itself. 

Comparing VEFI and LAPI orbital morphology was similarly unfruitful. The 
zeroth order expectations of both the VEFI and LAPI data showing two cuts through a high 




latitude auroral region per hemispheric pass is present But for each auroral crossing the 
structure in the precipitation and convection show little systematic variation or correlation. 
Comparing with the images leads to further complications. The LAPI fine structure in each 
auroral crossing is typically contained within 4 to 6 pixels which were observed within ± 
12 minutes of the LAPI data. This leads to a temporal and spatial mismatch that is not 
overcome. Basically, the variability is on time scales less than 12 minutes. 




2. CONCLUSIONS AND RECOMMENDATIONS 


The final objective of tins study has been found to be unachievable with present day 
data for the following reasons: 

• Present day auroral images do not have an adequate spatial nor temporal 
resolution. 

• Electric field data acquired every orbit (90 - 100 minutes) are uncorrelated. 
Geomagnetic and solar wind variability is on time scales much shorter. 

• Similarly, auroral precipitation is only weakly correlated from one pass to the 
next. 

Consequently, both the observed electric field and precipitation is changing in what 
is mathematically a step wise sense which precludes determining either the time constants 
of the change or specifically which geomagnetic, solar wind, or other indices are driving or 
correlated with them. Although this conclusion may have been a reasonable guess prior to 
the study it was also expected that the combined DE data base being the most 
comprehensive to date might shed new light on the problem. Based upon this study, 
specific requirements can be imposed upon future data sets such that an improvement in our 
knowledge of the large scale magnetospheric electrodynamics can be made: 

• Images need to have both enhanced spatial ( ~ 10 km) and temporal resolution 
(10’s of seconds). The POLAR satellite will produce such a data base. 

• These images need to be complete sequences spanning several DE-2 type orbit 
periods (hours). 

• In situ particle and field measurements do not contribute significantly to this 
problem unless they can be made at multi-locations with repetition rates on the 
order of minutes. 

• The (fine) structure observed in a single present-day pass does not contain 
information which can be correlated to activity indices or other indices. This is 
readily proven when a single pass is compared to a statistical model output for 
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the appropriate activity indices. The agreement is poor, the standard deviations 
and mean values are comparable. 

• There is a need for more electric field-current measurements being acquired at 
the same time with repetition rates on the same time scales as images (10’s of 
seconds). The Assimilate Mapping of Ionosphere Electrodynamics (AMIE) is a 
technique which combines various electrodynamic data sets into a simple global 
electric field and current distribution. However, to date, this technique is 
heavily data starved making its predictions statistically questionable. 

• The data to answer this problem requires the deployment of many satellites. A 
constellation of 50 small electrodynamic satellites has been put forward by 
Sojka [1989] and Sojka and Redd [1990]. 

This requirement for many satellites described by Sojka [1989] and Sojka and Redd 
[1990] is a direct result of the research earned out under this grant 
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Table 1. Geomagnetic conditions 


Day Number 

Date 

Average Kp 

Designation 

312 

November 8, 1981 

4 

active 

326 

November 22, 1981 

2+ 

moderately active 

327 

November 23, 1981 

3+ 

active 

331 

November 27, 1981 

1 

quiet 

333 

November 29, 1981 

1 

quiet 

351 

December 17, 1981 

1+ 

quiet 

356 

December 22, 1981 

1- 

quiet 

365 

December 31, 1981 

3 

moderately active 
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Table 2. UV image data base 


Day Number 

Number of Image Passes 

Total Number of UV Images 

312 

3 

41 

326 

2 

28 

327 

3 

50 

331 

2 

33 

333 

1 

17 

351 

3 

54 

356 

2 

41 


365 


3 


49 
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L Introduction 

The author has taken a parochial view of the subject matter 
under the title of ionospheric physics. Physics is used as the 
key to separating this section of the aeronomy review from the 
other aeronomy reports. Specifically, work categorized as 
chemistry, optical emissions, and thermosphere will not be 
reviewed in this section. However, coupling to the 
thermosphere and more specifically to the magnetosphere is 
considered. A brief section is also included on active 
experiments which is a productive area in plasma physics and 
also to some extent in ionospheric physics, but one that has 
fallen on hard times with regard to the U.S. funding agencies. 
To many, the ionosphere and its physics is regarded as a 
mature applications science or perhaps even a technology. 
This is far from the truth. The original wave of zeroth order 
exploratory experiments and theories have long since gone, 
but we still do not have predictive models (theories) that are 
worthy of being called technological tools. There is 
considerable current interest in developing a capability to 
predict the morphology of ionospheric parameters on a global 
scale. Some shortcomings for achieving this goal are 
summarized in Section 4, the conclusion; they define, in part, 
the future research goals of the ionospheric physics 
community. Section 2 divides ionospheric physics according 
to the community's normal regional classification, i.e. 
equatorial, mid-latitude, auroral, and polar ionosphere. In 
addition, the mid-latitude is extended upwards to include the 
cold plasmaspheric population, and the polar region is also 
extended upwards to include the polar wind. To date, this polar 
wind population has still not been fully surveyed. Not all 
ionospheric physics is naturally restricted to a regional 
discussion and, therefore, in Section 3 the topics of 
ionospheric electrodynamic inputs, ionospheric irregularities 
and waves, active experiments, forecast modeling, and 
coupling to other regions are reviewed, but material already 
discussed in Section 2 is not repeated. 

A measure of the balance of research in ionospheric 
physics is obtained from the distribution of experimental, 
analytical, and theoretical publications that have appeared 
during this reporting period. Of course, the author has applied 
subjective criteria in classifying papers according to these 
three categories. However, the results are still representative 
of the balance in the U.S. ionospheric community. Figure 1 
shows this distribution not only by research category but also 
by the four regional areas identified in Section 2. The data used 
to create Figure 1 are taken from the Journal of Geophysics 
Research and Geophysics Research Letters in the reporting 
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period up to August, 1990, including a total of 172 papers by 
U.S. authors (at least one U.S. author was the criterion for 
inclusion in this review). As expected, auroral physics 
dominates the research effort because it includes the most 
diverse collection of dynamic and exotic effects, and it is also 
in this region that coupling with both the thermosphere and 
magnetosphere is most prominent. Currently, the polar 
ionosphere is receiving renewed attention especially in its 
coupling to the open magnetosphere/solar wind for northward 
IMF* conditions. Several ionospheric physics review articles 
have appeared since the last IUGG report; papers by Heelis 
[1987], Killeen [1987], and Inan [1987] are the most relevant 
to this review. Fraser-Smith [1987] reviewed the magnetic 
field topology in the ionosphere while Gorney [1990] 
discussed the solar cycle effects on the near earth 
environment. Measuring wavelength in space plasmas was 
considered by LaBelle and Kintner [1989] and high latitude F 
region plasma irregularities were reviewed by Tsunoda [1988]. 
LaBelle [1989] reviewed studies of auroral radio noise during 
the past twenty years. A series of review articles described the 
contribution of the first five years of the Dynamics Explorer 
(DE) program, of which the most relevant to this review are 
those by Hoffman [1988], Heelis [1988], and Burch [1988]. 
Demars and Schunk [1987] reviewed temperature anisotropies 
in the terrestrial ionosphere, Schunk [1988] reviewed the 
theory and numerical implementation of ionospheric models, 
and Sojka [1989] reviewed the status of global scale, physical 
models of the ionosphere. 

2.1. Equatorial Ionosphere 

Research in this region focused on the observation and the 
analysis of plasma density, composition and drift, equatorial 
bubbles, and instabilities. Further analysis of the CONDOR 
experiments on E region instabilities was carried out. 
Modeling and theory focused on large scale equatorial plasma 
and instability studies. The equatorial F-layer, especially in 
the night sector showed large scale structure, i.e., the 
equatorial anomaly. Our understanding of these structures is 
qualitatively good but quantitatively poor. These structures 
arise from a combination of E x B and neutral wind induced 
drifts. Much of the work focuses on separating these drift 
effects and quantifying them. Eastward electric field data from 
Dynamics Explorer (DE)-2 (Aggson et al., 1987, Maynard et 
al., 1988, and Coley and Heelis, 1989) together with 
incoherent scatter radar drift data from Jicamarca (Pingree and 
Fejer, 1987 and Fejer et al., 1989) were used to observe the 
equatorial plasma drifts. These analyses show a high degree of 
compatability and allude to the strong local structure present, 
especially in the evening post sunset region. Fejer et al., 
[1990] showed that during disturbed storm periods equatorward 
penetration of high latitude electric fields is preferentially in 
the zonal component. Night enhancements of 630 nm 
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Fig. 1. Relative distribution of ionospheric physics publications for JGR and GRL in this reporting period as a 
function of research methodology and ionospheric latitudinal region. 


emission observed by the AE-E satellite led Fesen and Abreu 
[1987] to conclude that the local time distribution of plasma E 
x B drift needed to be offset by an hour from the normally 
adopted convection pattern. During the main phase of strong 
storms Sahai et al., [1988] observed 391 and 486 nm 
emissions and inferred the presence of energetic particle 
precipitation and used the 630 nm emissions to follow the F- 
layer vertical motions during storm periods. Campbell and 
Schiffmacher [1988] analyzed magnetometer data for 1965 to 
obtain hemispheric asymmetries, seasonal, and kp dependence 
of the Sq current system. 

Bubbles in the equatorial plasma were studied by rocket 
experiments which were able to create them chemically 
[Klobuchar and Abdu, 1989] and Rohrbaugh et al., [1989] were 
able to track them optically along magnetic field lines from 
Hawaii. The instabilities present in the equatorial plasma were 
remotely sensed by radar techniques: a 50 MHz radar 


interferometer [Kudeki and Farley, 1989], a 250 MHz signal 
from geosynchronous satellites [Bhattacharyya et al., 1989], 
and a combination of VHF and HF radars [Kudeki et al., 1987]. 
The latter research was in conjunction with the CONDOR 
project in which rocket instruments were flown through the 
equatorial E region. In situ data from this project led to direct 
observation of the plasma structure associated with plasma 
instabilities [Pfaff et al., 1987 a,b]. These papers focused on 
the analysis of gradient drift and two stream waves. Ronchi et 
al., [1989] carried out a linear global eigenmode analysis of 
the gradient drift instability. The effect of small scale 
turbulence ( clOOm) on irregularity models of the gradient 
drift instability of kilometric scale was developed by Ronchi 
et al., [1990]. Maruyama [1988] developed a diagnostic model 
of linear growth rates for the gravitational Rayleigh-Taylor 
instability. These models allow for the effects of E x B and 
wind induced drifts and, hence, have application to equatorial 
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spread F and equatorial bubble phenomena. The generation of 
gravity waves by a solar eclipse was observed using TEC 
measurements and was found to be limited to observing 
stations within 500 km of the eclipse path [Singh et aL, 
1989]. Klobuchar and Lee [1989] observed amplitude 
variations in VHF signals as a precursor to onset of strong 
equatorial scintillations, although they found no modulation 
of the TEC measurements. 

Satellite composition measurements were used to study the 
equatorial plasma. Grebowsky and Reese [1989] used AE-E ion 
spectrometer measurements of Fe+ to study the diurnal and 
seasonal equatorial F region distribution. Their analysis 
found trends which were not entirely consistent with earlier 
theoretical studies. Heelis et al., [1990] used DE-2 plasma data 
to study the He + distribution with special emphasis on 
seasonal trends and a comparison with mid-latitude data. The 
response of the equatorial ionosphere to neutral atmosphere 
perturbations induced by magnetic storms was modeled by 
Fesen et al., [1989]. Separate ionosphere and neutral 
atmosphere (TGCM) modeling was carried out to simulate the 
March 22, 1979 storm; agreement between plasma 
calculations and observations was limited to the 140° W 
sector. Sojka and Schunk [1989] carried out a study of the 
seasonal dependence of the global ionosphere. They found 
significant seasonal (and hemispherical) differences in the 
Appleton anomaly and attributed these to asymmetries in the 
neutral atmosphere. A fully analytic, low and middle latitude 
ionosphere model was developed by Anderson et al., [1989]. 
This model has seasonal a_nd solar cycle dependence and 
reproduces results of earlier semi-empirical models. 

2.2. Mid-Latitude Ionosphere And Plasmasphere 

A wide range of experimental, theoretical, and analytical 
research techniques has been used in mid-latitude ionospheric 
physics. Although these techniques are not new, 
collaboration in coordinated campaigns has resulted in 
significant progress being made in this area. The upward 
extension of the mid-latitude ionosphere into the 
plasmasphere as well as downward coupling from the 
magnetosphere in the form of stable auroral red (SAR) arcs has 
similarly benefited from this coordinated research approach. 
The extension of the F region mid-latitude trough into sunlight 
was extensively studied by Whalen [1987, 1989]. He showed 
that this feature was spatially well-defined and dependent, both 
on longitude and magnetic activity. Sojka et al., [1990] 
showed that these empirical relationships were present in their 
time dependent ionospheric model. Heelis et al., [1990] using 
DE-2 plasma composition data, showed that at 900 km He + is 
the dominant ion in the mid-latitude night sector and argued 
why this He + dominance was restricted to mid-latitudes. The 
presence of Ca + , its altitude distribution, and relationship to a 
meteor shower [Ursids, (December)] was studied by Torr et al., 
[1990]. They observed a resonance fluorescence emission 
from Ca + , inverted the surface brightness profile, and produced 
volume emission rate profiles. Emission peaks were present 
in both the E and F region during the shower. Rasmussen et 
al., [1988] compared ionospheric densities and electron 
temperatures from incoherent scatter radar observations 
(Chatanika, Millstone Hill) from a MITHRAS study with the 
TDIM. They found that although the densities throughout the 
F region were in agreement, the electron temperature 


interpretation was ambiguous. The mid-latitude F region 
density and height were studied using data from 41 iono sonde 
stations from a SUNDIAL campaign and fitted with TDIM 
simulations by adjusting the neutral wind [Sica et al., 1990]. 
This study showed that the topside 0+ flux was a significant 
unknown parameter, preventing a unique wind solution from 
being obtained. The solar cycle dependence of the neutral wind 
at both Boulder and Wallops Island, USA, was deduced from 
ionosonde data by Buonsanto [1990]. Sojka et al., [1988] 
compared the "quiet day" diurnal, seasonal, and solar cycle 
variations of Fof2 at Argentine Islands with the TDIM and 
obtained very good agreement, Berkey and Stonehocker 
[1989] compared different methods of true height analysis of 
ionosonde data. The effect of conjugate ionospheric 
photoclectrons upon the net ionization rates of 0 + and N2 + 
was modeled as a function of altitude by Richards and Torr 
[1988]. They found that these ionization rates were a function 
of altitude which is complicated by the presence of the 
conjugate photoelectrons down to around 250 km. Sojka and 
Schunk [1989] using the TDIM showed that the hemispherical 
difference in the mid-latitude seasonal anomaly was associated 
mainly with the hemispheric difference in the neutral 
atmosphere. Detailed longitudinal dependencies and absolute 
ratios of the seasonal anomaly were, however, only 
qualitatively modeled. Anderson et al., [1989] developed a 
fully analytic low mid-latitude ionospheric model. 

The dynamics of the mid-latitude ionosphere, especially 
during storm periods, is driven by enhanced electric fields, 
neutral drifts, and composition changes. Fejer et al., [1990] 
using incoherent scatter radar data from the January, 1984 
GISMOS campaign, studied the penetration of high latitude 
electric fields through the mid-latitude region and showed that 
there were differences between the observation and the theory. 
Burnside et al., [1987] improved the Arecibo plasma drift 
measuring technique and Berkey et al., [1990] analyzed the 
solar cycle dependence of the plasma drifts at Arecibo and 
found the most marked Nmf2 variability in the night sector. 
Maynard et al., [1988] using DE-2 electric field observations, 
computed average meridional electric fields at mid-latitudes. 
Although various dependencies were found, no consistent 
magnetic activity trend was found in these averaged data. The 
effects of major storms on the ionosphere were studied in 
several different ways; Forbes et al., [1988] used a harmonic 
analysis of F region data from a latitudinal chain of 
ionosondes to show ionospheric storm effects; Hajkowicz 
and Hunsucker [1987] used the periodic variation of the F- 
layer's virtual height in both the northern and southern 
hemisphere to locate the source of TID's at the foot of the 
energetic particle precipitation belts; Richards et al., [1989] 
using Millstone Hill ETS data, found that vibrationally excited 
N2 did not account for the ionospheric negative storm phase 
and that during the storms, MSIS 86 densities were not 
reliable; while Yeh and Foster [1990] observed exceptionally 
large O+ outflows in the mid-latitude ionosphere above 
Millstone Hill during a large storm. Coherent-array HF 
Doppler techniques were developed to monitor TID's [Adams et 
al., 1988; Jacobson and Carlos, 1989 a,bj. Storm effects also 
manifest themselves in the thermosphere. The use of 
ionospheric observations to deduce neutral atmospheric 
parameters at Arecibo, was studied by Burnside et al., [1988] 
who showed the difficulty in obtaining agreement with MSIS 
[Richards et al., 1989 and Oliver, 1990]. From Millstone Hill, 
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Foster and Aarons [1988] observed strong, 300m/s, westward 
convection in the morning sector as a region of storm 
enhanced sunward convection retreated poleward. These 
observed westward drifts were equatorward of the 
magnetospheric shielding layer and were driven by the storm 
enhanced neutral winds modeled by Spiro et al., [1988], Hagan 
[1988] studied the effect of the February, 1986 storm on the 
thermosphere while Forbes and Roble [1990] investigated the 
F region storm effects on the thermosphere (TGCM) as a 
precursor to a fully self-consistent ionosphere-thermosphere 
study. Williams et al., [1987] discussed coupling of the D and 
E region with the middle atmosphere. The D and E region 
incoherent scattering observations were compared with theory 
by Ganguly and Coco [1987]. 

Rasmussen et al., [1988] developed a combined 
conductivity and E region bottomside F-layer (85-220 km) 
model and showed their E region model densities to be within 
50% of those observed at Chatanika and Arecibo incoherent 
scatter radar observations. They also found that at night the Fi 
densities can contribute as much as 80% to the height 
integrated Pedersen conductivity. Cole [1990] studied the 
electric currents flowing in the ionospheric E region and 
developed new expressions for the components of heat 
dissipation and conductivity associated with the ionospheric 
current drivers. He also derived a new diffusion equation for 
magnetic induction in E region like plasma regions. 
Buonsanto [1990] developed an E-Fi daytime photochemical 
equilibrium model and compared it with extensive mid-latitude 
data sets. He concluded that to get agreement with 
observations, the EUV fluxes needed to be increased by 25- 
30% and that at equinox MSIS 86 neutral densities needed to be 
decreased by about 25%. In winter, this neutral density 
decrease would have to be even larger. These shortcomings in 
the MSIS model were consistent with the Arecibo analysis of 
Burnside et al., [1988]. St-Maurice et al., [1989] analyzed 
coherent echoes at 44 MHz from the Millstone Hill radar. 
They observed these echoes in the 105 to 115 km region north 
of the radar where the viewing geometry was suitable for 
coherent scatter. At Arecibo, Djuth and Gonzales [1988] used 
the 430 MHz radar during high-power, high frequency 
experiments to study the growth of Langmuir waves. The 
conversion of HF into Langmuir waves was put forward by the 
authors as a promising mechanism for the production of 
Langmuir waves in sporadic E. Riggin and Kadish [1989] 
modeled the non-local growth rate of gradient drift plasma 
waves and compared those with local growth rates of kilometer- 
scale waves observed in the vicinity of mid-latitude sporadic E 
layers. However, they found that the non-local effects alone 
do not explain the dominance of the kilometer-scale waves. 
Fukao et al., [1988] observed lm small-scale field-aligned 
irregularities in the mid-latitude F region using the MU (VHP) 
radar. Their results showed that the irregularities were, on 
occasions, quasi-period on both space (50 km) and time (12 
minutes) scales. Rodger and Aarons [1988] compared 
conjugate occurrence frequency of irregularities and deduced 
both seasonal and storm correlations. Djuth et al., [1987] 
observed enhanced ion-acoustic waves (lm wavelengths) with 
a 138.8 MHz backscatter radar at Arecibo during the high- 
power, high frequency F region modification experiments. 
Hansen and Henley [1988] modeled the effect of isotropic 
turbulence in a background F region plasma into which a 
barium release was made. They found that the turbulence 


enhanced the diffusion coefficient for the release and, hence, 
could explain the predominance of large length scale 
turbulence observed during such releases. 

Slater et al., [1987] found enhanced 10 eV electron 
precipitation (DE-2 data) in conjunction with 23 instances of 
SAR arcs observed by ground based photometers. The electron 
energy fluxes were found to be sufficient to produce the 630 nm 
SAR emission rates. Kozyra et al., [1987] using both DE-1 
and DE-2 correlated ring current heavy ion fluxes at 17 keV 
with SAR arc electron temperature enhancements. Such 
energetic heavy ions (0 + ) provide the electron heating (10 eV 
electrons) needed for the SAR arc through Coulomb collisions. 
Mendillo et al., [1989] studied three distinct sub-auroral 630 
nm emission features, detached arcs, ripples in the diffuse 
aurora, and patches. DMSP satellite and Millstone Hill radar 
data were used to infer the sources of these emissions. The 
ripples in the diffuse aurora appear to be driven by enhanced 
ion heating due to enhanced plasma convection while the arcs 
(SAR) were associated with F region density enhancements and 
ion precipitation, and the patches appear to be associated only 
with enhanced radar returns from the E region. Okano and Kim 
[1987] used photometric triangulation to infer the SAR arc 
emission height to be 370-460 km on the 26/27 September, 
1979. This height rapidly decreased as the SAR arc brightened 
in response to substorms. The relative effects of the 
ionosphere-thermosphere and the magnetosphere on SAR arc 
formation was studied by Kozyra et al., [1990] as a function of 
season and solar cycle. The seasonal asymmetries, especially 
around solstice, can best be explained by the ionospheric- 
thermospheric asymmetries. However, the solar cycle 
dependence can best be explained by a reduction in 
magnetospheric heat flux from solar maximum to solar 
minimum which is attributed to the compositional change in 
the magnetospheric plasma. 

The strong interdependence of the plasmasphere and the 
mid-latitude ionosphere justifies regarding the plasmasphere 
in the context of ionospheric physics. Brace et al., [1988] 
used both DE-2 F region and DE-1 plasmasphere electron 
temperatures to correlate latitudinal boundaries in these 
regions. In the 2200 to 1300 LT sector (i.e., midnight, 
morning, and noon), the Te signatures are essentially 
identical. Outside this sector, the plasmasphere bulge 
complicates the correlation since no similar feature is found in 
the ionosphere. In the region of electron heating at the outer 
plasmasphere, Roberts et al., [1987] found enhanced heavy 
ion [0 + , 0 ++ , N + ] densities. The latitude of the heavy ion 
enhancements tracked the magnetic activity variation (Dst) 
and was associated with an ionospheric heavy ion outflow. 
Topside H + fluxes were not found to follow the same diurnal 
variation as the 0 + fluxes at Arecibo, but did show flows 
during magnetic storm disturbances [Tepley and Ken, 1989]. 
At higher altitudes in the plasmasphere, Menietti et al., [1988] 
studied the ion flow morphology using DE-1 data. Horwitz et 
al., [1990] used DE-1 and 2 composition data and the FLIP 
model to compare the composition and density of the F region 
and the plasmasphere. The 0+ latitudinal variation in the 
plasmasphere correlated with that of 0 + in the F region, but 
not with the light ion H+. The He+/H+ ratio was greater than 
unity in the topside ionosphere, but was around 0.2 in the 
plasmasphere. When the FLIP model allowed for fractional 
trapping of ionospheric photoelectrons, the results of 
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calculations were in good agreement with the observations. 
Horwitz et al., [1990] used the RIMS database to characterize 
both the plasmasphere's density structure and boundary 
locations. Singh and Torr [1990] used a hydrodynamic model 
to study the effect of ion temperature anisotropies on 
interhemispheric plasma transport. Rasmussen and Schunk 
[1990] developed a three-dimensional time dependent, non 
-linear hydrodynamic model of the plasmasphere which 
included cross-L-shell drifts. The initial solar minimum study 
of equatorial densities gave results within a factor of 2 of the 
average whistler and satellite density observations. The mid 
-latitude ionosphere is at the feet of the closed magnetic flux 
tubes which duct VLF signals, and this same magnetic field 
geometry supports the ring current wave-induced precipitation 
events; both the ducting of the VLF signals and the wave 
-induced precipitation were the topic of several studies 
[Neubert et al., 1987; Carpenter et al., 1988; Carpenter and 
Orville, 1989; Burgess and Inan, 1990; Jarvis et al., 1990]. 

2.3 . Auroral Ionosphere 

The auroral ionosphere contains a diverse collection of 
phenomena; very few of these can be understood in isolation. 
No clear distinction between magnetospheric and ionospheric 
phenomena exists here. Dynamics are present both spatially 
and temporally over a wide range of scales. Work in this field 
continues to be prolific and advances often appear to further 
complicate rather than resolve problems. In the following 
paragraphs the research is reviewed by forcing a subjective 
ordering, but subsections are not defined as this would be a 
futile categorization. 

Both theoretical and empirical research has found the 
auroral ionosphere abound in small-scale irregularities and 
waves. The classification and subsequent understanding of 
instabilities has led to a continuing series of studies. In turn, 
our knowledge of the plasma structuring mechanisms has 
advanced significantly. Observations in the F region by 
Ruohoniemi et al., [1987] using HF coherent scatter and 
incoherent scatter radars from Goose Bay and Sondrestrom 
Fjord, respectively, showed that small-scale (13.9m) field 
aligned irregularities drift at the ambient plasma drift speed. 
Using the same HF radar, Ruohoniemi et al., [1988] observed 
similar coherent backscatter from decameter scale plasma 
density irregularities for periods of about 1 hour in the evening 
sector at sub-auroral latitudes. Foster et al., [1988] found that 
certain Millstone Hill radar spectra which were typically 
regarded as n hard-target" backscatter were also consistent with 
intense radar returns associated with topside unstable plasma 
regions of current driven ion acoustic modes. Using rocket 
data, Earle et al., [1989] found that in the presence of large 
velocity shears in the F region there were both turbulence and 
electrostatic waves and concluded that the velocity shears were 
associated with field aligned currents which drove the 
electrostatic ion cyclotron waves. Theoretical work on the F 
region also covered several topics. Chaturvedi and Huba 
[1987] looked at the Interchange instability associated with 
plasma blobs (300-600km) and found that it could generate the 
scintillation causing irregularities (1.1 km wavelength) 
associated with the blobs. Cornwall [1988] found exact 
solutions for the Liouville equation of an adiabatic auroral arc. 
These solutions gave stability limits for processes such as the 
Kelvin-Helmholtz instability. Basu and Coppi [1988, 1989] 


developed a theory of collisional electrostatic modes driven 
unstable by sheared field aligned ion flow velocities and 
showed that this theory could predict the observed [Basu et al., 
1988] fluctuations from several meters to a kilometer. 
However, Ganguli et al., [1989] argued that the transverse 
velocity shear would readily stabilize the above unstable 
modes. Nishikawa et al., [1990] modeled electrostatic 
turbulence due to both the transverse and parallel sheared 
flows. The effect of ion-ion and ion-neutral collisions on 
quasi-linear ion heating by the current- driven ion cyclotron 
instability was studied by Satyanarayana et al., [1989]. In the 
lower F region (200-300km), they found that the collisions 
isotropize the temperature distribution while at higher 
altitudes (600km) anisotropies can be sustained even with the 
collisions. 

Experimentally, Noble et al., [1987] studied E region 
irregularities of 1 to 7m using the HF ionospheric heater 
facility near Tromso in conjunction with backscatter radars. 
Their results indicated that the E region irregularities were 
associated with instabilities at the upper hybrid resonance 
level rather than the reflection level. Valladares et al., [1988] 
used the Sondrestrom Fjord incoherent scatter radar to observe 
enhanced plasma lines from the topside of the particle 
produced E-Iayer. Villain et al., [1990] observed decameter 
irregularities of 10-20m jointly with the Goose Bay and 
Schefferville HF radars. They found evidence of two 
irregularity layers in the E region which are related to the 
action of the gradient drift and ion acoustic instability 
mechanisms. The Farley-Buneman instability was 
theoretically studied by Chaturvedi et al., [1987] and Machida 
and Goertz [1988]. Chaturvedi et al considered the effect of a 
parallel current on two-stream plasma instabilities and found 
that it modified the threshold for the onset of the Farley- 
Buneman instability, while Machida and Goertz studied the 
nonlinear saturation of the Farley-Buneman instability using a 
collisional 2.5 dimensional electrostatic particle simulation 
and found that the dominant mechanism for electron heating 
was an enhanced effective collision frequency rather than 
heating from the parallel electric field component of the 
waves. Kuo and Lee [1988] found that E region irregularities 
greater than 13m could be produced through a thermal 
instability causing filamentation of the auroral electrojet 
current. Farley and Providakes [1989] reconsidered the 
validity of the isothermal assumption in the theory of unstable 
ionospheric two-stream waves. They found that although the 
electrons are adiabatic, the ions are not. This leads to a 20 to 
40% difference in modeling the wave velocities with the newer 
calculations being in better agreement with the EISCAT and 
CUPRI observations. Satyanarayana et al., [1987] studied 
theoretically the effect of Coulomb collisions and parallel 
electrodynamics on the transverse Kelvin-Helmholtz 
instability, while Ganguli and Palmadesso [1988] studied the 
effect of parallel currents and transverse electric fields on 
electrostatic ion instabilities. 

Using the EISCAT incoherent scatter radar, Lockwood et 
al., [1987] observed bursts of intense poleward ion flows in 
the dayside auroral region. When the ion drift velocity 
exceeded the neutral thermal speed, they observed non- 
Maxwellian ion velocity distributions. Senior et al., [1987] 
observed both the E and F region density in the evening sector 
auroral oval using the Chatanika incoherent scatter radar data 
and found that the auroral E region extended 2° further 
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equatorward than the F region. Using both DE-2 and NOAA-6 
particle data, they showed that the extended E-layer was formed 
by proton precipitation in contrast to the higher latitude E-F 
layers which were formed mainly by electron precipitation. In 
the high latitude dawn sector, Robinson et al., [1988] found 
discrete arcs in diffuse auroral regions. The observed 
(Sonderstrom Fjord radar) electron density of the arcs peaked at 
180km while the diffuse auroral electron density peaked at 115 
km. Using both DMSP and HILAT satellites, they showed that 
the discrete arcs were associated with enhanced fluxes of low 
energy electrons. Kondo et al., [1990] statistically analyzed 
upflowing ion beam and conic distributions observed on DE-1 
and found them to be mostly associated with auroral field lines, 
although there was not a strong correlation with the direction 
of the field aligned currents. They also found that 0 + was 
accelerated more efficiently at higher levels of magnetic 
activity (Kp). On a global scale, Sojka and Schunk [1989] 
modeled the high latitude ionosphere using their TDIM model 
and showed the strong seasonal and UT control of the auroral 
ionosphere. Sojka et al., [1989] carried out a TDIM auroral 
substorm simulation using both the DE-1 and DE-2 data to 
define the auroral dynamics. The DE-1 images defined the 
global scale precipitation dynamics for the substorm. Their 
simulation indicated the need for a storm representation of the 
convection electric field. 

Schunk et al., [1987] used the TDIM model and the electron 
energy equation to study the effect of ionospheric return 
currents on the electron temperature. The uncertainty of the 
magnitude and direction of the magnetospheric heat flux was 
shown to be a major problem. Lilensten et al., [1990] used 
both Viking and EISCAT data to solve the electron energy 
equation along an auroral field line which enabled them to 
deduce the topside heat flux. Chiu et al., [1988] used a Monte 
Carlo simulation of wave particle ion heating in the auroral 
topside ionosphere to study the efficiency of the topside as a 
provider of magnetospheric plasma. They found it was very 
sensitive to charge exchange collisions and downward parallel 
electric potential drops. In addition, they predicted the 
presence of a heated (<10 eV) neutral population in the topside 
ionosphere. 

The auroral electron density profile is dependent upon the 
spectral shape of the incident precipitating particles, and the 
auroral emission is dependent upon the photochemistry of the 
auroral ionosphere. Jones et al., [1987] used electron density 
profiles measured by a radar to deduce the average energy of the 
incoming electrons and then compared this to the 1(630.0) 
/I(427.8) and 1(844. 6)/I(427.8) optical emission ratios (see 
also Niciejewski et al., [1989]). Rees et al., [1988] deduced 
auroral energy deposition rates, characteristic electron energy, 
and ionospheric parameters from DE-1 auroral images at three 
wavelengths. The ionospheric conductance was deduced from 
the DE-1 auroral images and was checked against the 
conductance inferred from Chatanika incoherent scatter radar 
data [Robinson et al., 1989]. All sky television at Syowa 
station and particle data from the DMSP-F6 satellite were used 
by Ono et al., [1989] to analyze the dusk auroral emissions 
caused by protons and electrons. Theoretical work on the ratio 
of auroral emissions to obtain the characteristic energy of the 
electron spectrum was done by Rees and Lummerzheim [1989] 
who found the I(337.1)/I(427.8) ratio to be better than the 
I(630.0)/I(427.8) ratio. Richards and Torr [1990] modeled the 
337. lnm emission rate dependence upon the characteristic 
energy of the auroral electrons. Morrill and Benesch [1990] 


studied how ’’preconditioning" the N2 gas affects the 
production of excited N2 vibrational distributions and 
subsequently, the N2 band intensities observed in the aurora. 
Enhanced soft auroral precipitation (clkeV) was observed by 
rocket borne instrumentation in the morning sector auroral 
zone in conjunction with electron density enhancements 
[LaBelle et al., 1989]. These enhanced densities were modeled 
and were found to be produced by the soft precipitation. Basu 
et al., [1987] compared a linear transport theory model of 
auroral ion (proton) precipitation with Chatanika radar 
electron densities and the NOAA-6 proton precipitation data 
and found good agreement. Hardy et al., [1989] binned the 
auroral ion precipitation measured by DMSP-F6 and F7 
satellites as a function of Kp to produce a statistical ion 
precipitation model. The dependence of ionospheric 
conductivity calculations upon the flux and the energy of the 
precipitating electron spectrum was used by Robinson et al., 
[1987] to highlight the difference between the characteristic 
energy of a Maxwellian distribution versus the average energy 
of an auroral spectrum. Fuller-Rowell and Evans [1987] 
produced a model of the height integrated Pedersen and Hall 
conductivities from TIROS-NOAA satellite data, while Hardy et 
al., [1987] used DMSP data to model the pattern of ionospheric 
conductivity as well as auroral energy and number flux. Foster 
et al., [1989] combined the Fuller-Rowell and Evans [1987] 
conductivity model with Millstone Hill convection patterns to 
deduce auroral ionospheric currents. Robinson et al., [1989] 
developed a procedure for the inversion of auroral 
Bremsstrahlung x-ray spectra to obtain electron deposition 
information, while Ahn et al., [1989] used such a technique to 
deduce auroral conductivities which they combined with ground 
based magnetometer data to infer global scale heating rates. 
Kamide et al., [1989] used ground magnetometer data together 
with DE-1 and 2 observations to deduce ionospheric 
electromagnetic quantities. 

Jet stream wind models were shown to produce electric 
fields parallel to the magnetic field with current closure in the 
F region [Nakada, 1987]. Additionally, Nakada [1989], using 
an oscillatory electrical circuit, showed that the jet stream 
could cause oscillations in the Pc-1 to Pc-3 range which is in 
accordance with observations. Poulsen et al., [1990] 
developed a three-dimensional night time D region model to 
study the effect of localized perturbations upon sub- 
ionospheric VLF propagation. 

The cleft/cusp region has received attention due to its key 
location in relation to plasma flow into the polar cap as well 
as mapping FTE signatures. Sandholt et al., [1989] used 
HILAT plasma data and ground observations of polar cusp 
auroras to deduce the electrodynamics of the cusp. Using 
ground-based search coil magnetometer, riometer, photometer 
and ELF-VLF receiver data from South Pole and McMurdo 
stations, Engebretson et al., [1990] looked for Pc-3 pulsation 
signatures in the ionosphere. They found low energy 
precipitating magnetosheath-like electrons in the cleft which 
were, at times, modulated with frequencies similar to those of 
upstream waves. Wolfe et al., [1990] used conjugate stations 
at Igaluit and South Pole to study Pc-3 pulsations near the 
cusp. Very few conjugate and simultaneous observations of Pc- 
3 were made. The authors concluded that upstream ion 
cyclotron waves were the source of the Pc-3 events. Greenwald 
et al., [1990] observed simultaneous conjugate dayside 
convection and noted rapid response to changing IMF-By 
conditions. Yeh et al., [1990], using the Millstone Hill 
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incoherent scatter radar, observed the dayside cusp extending 
as low as 60° invariant latitude during the large February 8-9, 
1986 storm. 

The FTE electrodynamic signatures are associated with the 
cusp region. Lockwood et al., [1988] and Smith and 
Lockwood [1990] used the EISCAT radar to observe the plasma 
flow in the cusp and inferred convection signatures consistent 
with a simple twin-vortex FTE model. Strong upward plasma 
flows were observed by AMPTE-UKS in these events [Smith 
and Lockwood, 1990]. Saflekos et al., [1990] indirectly 
observed evidence of reconnected flux tubes which they 
interpreted to be extensions of FTE’s found near the 
magnetopause. The FTE signatures were studied theoretically 
by Zhu and Kan [1989]; the ground magnetic signature was 
studied by Wei and Lee [1990]; the ground magnetic signature 
for a moving elongated plasma cloud was modeled by Crooker 
and Siscoe [1990]; and the convection footprint in the 
ionosphere was studied by Crooker [1990]. It is in the 
ionosphere that high resolution two-dimensional analysis of 
the FTE can be made, and several studies have focused on the 
expected and observed FTE ionospheric signature [McHenry 
and Clauer, 1987; Bering et aL, 1988; Kan 1988; Basinska et 
al., 1989; Lockwood and Smith 1989]. 

Briining et al., [1990] and Potemra et al., [1990] used 
Viking satellite data to study dynamic discrete auroral 
structures in the 1400 MLT time sector. Bruning et al., 
considered the current systems and parallel potential drops and 
found evidence for a saturation current in intense postnoon 
aurora, while Potemra et al., studied the temporal variation of 
"bead” or "pearl" auroral forms from the UV images. Using 
particle distributions, waves, and dc electric and magnetic 
fields measured from two sounding rockets, McFadden et al., 
[1990] determined the electrodynamics of a thin (10 km) 
auroral arc. How the neutralization of the auroral arc upward 
current depends on arc thickness was theoretically studied by 
Nakamura et al., [1989] and how the upward currents are 
distributed across the arc was theoretically studied by 
Greenspan [1989]. Porter et al., [1987] developed a 
computationally fast electron transport deposition code using 
a Neuman iteration method. Reiff et al., [1988] used 
simultaneous high altitude (DE-1) and low altitude (DE-2) 
satellite data to observe indirectly the auroral parallel 
potential drop, while Kan and Cao [1988] considered the effect 
of an auroral parallel potential drop on magnetosphere- 
ionosphere coupling. Observations were also made to extend 
our empirical understanding of the ionospheric role in 
generating auroral signatures [Evans et al., 1987; Rosenberg 
et al., 1987; Imhof et al., 1988; Miah 1989]. Sandahl et al., 
[1990] studied the correlation of auroral precipitation (Viking 
UV images) and convection (EISCAT) in the midnight auroral 
region. During the substorm breakup when the auroral 
poleward boundary expands poleward, Robinson et al., [1990] 
found that this motion did not correlate with the variation in 
the local convection electric fields. Boehm et al., [1990] used 
high resolution sounding rocket observations to infer the 
presence of shear Alfven waves with amplitudes >100 mV/m 
in the auroral zone. The auroral zone plasma convection was 
modeled empirically by Holt et al., [1987]. They used 
Millstone Hill radar observations and compared their model 
with earlier results; however, in the average model discrete 
features like the Harang discontinuity are lost. Sica et al., 
[1988] modeled theoretically the high latitude ionosphere to 
find conditions under which large Field aligned ion drifts could 


be maintained. 

Temporal storm and substorm features associated with the 
auroral zone ionosphere were also studied. Gledhill et al., 
[1987] used data from Sanae, Antarctica to study the 
electrodynamics of the substorm and found that the 
magnetometer D component gave a positive spike instead of 
the normal negative spike. Robinson and Vondrak [1990] 
deduced the electrodynamic properties of auroral surges using 
Chatanika and DE-1 and 2 data. Similarly, Lyons et al., 
[1990], used Viking satellite, magnetometer, and radar data to 
study auroral surges in the vicinity of the Harang 
discontinuity. Pangia et al., [1990] correlated Pc-5 magnetic 
pulsations with substorm onset using Goes-2 and 3 
geostationary data in conjunction with ground magnetometer 
data. The AMIE technique was used by Knipp et al., [1989] to 
deduce the electrodynamic patterns during the ETS disturbed 
day of September 19, 1984. Strong eastward convection in the 
midnight sector of the auroral oval during the first large 
substorm was found to be associated with a marked westward 
rotation of the two-cell convection pattern. Zhu and Kan 
[1990] considered theoretically the effect of ionospheric 
recombination time scales on auroral substorm signatures. 
Their model produced discrete aurora along the poleward 
boundary of the diffuse aurora in the evening-midnight sector. 
Stern [1990] analyzed qualitatively the auroral substorm using 
a one-dimensional model. The critical dependence of Region II 
field aligned currents upon ionospheric conductivity gradients 
was highlighted by Kan [1987] and Cheng et al., [1987]. 

Finally, a few auroral papers that do not fit into the above 
topic catagorization by paragraphs. Hones et al., [1989] used 
DE-1 auroral images and found that for about one third of the 
time, when IMF Bz was northwards, the auroral oval was 
configured like a "horse collar". This pattern often evolved 
into a theta-like auroral pattern. Robinson et al., [1989] 
studied the plasma and field properties of auroral suprathermal 
electron bursts and found that the bursts, typically < IkeV were 
strongly correlated with enhancements in the ac electric field 
at frequencies below 1kHz. This led the authors to suggest that 
the electrons were heated via interaction with ion cyclotron 
waves. Hultqvist et al., [1988] using Viking data, observed 
both energetic electrons and ions of similar energy moving 
upward along magnetic field lines and concluded that these 
Field aligned fluxes were associated with strong broadband low- 
frequency turbulence and a downward current. They inferred 
that the waves accelerated the electrons (a resonance effect), 
but that the ions "see” it only as a wave field (no significant 
heat extraction). Calvert and Hashimoto [1990] studied 
magnetoionic modes and the propagation properties of auroral 
radio emissions using the DE-1 radio spectrograms and then 
ray traced their propagation behavior. They found that all four 
modes which could propagate were present in the data and all 
were associated with specific regions in the auroral zone. Sa 
[1990] observed VLF emissions triggered by discrete waves 
transmitted from Siple Station and received at Lake Mistissini, 
Quebec. Killeen et al., [1988] used DE-1 and 2 data of the 
auroral and polar ionosphere to study the relationship between 
the auroral morphology and the high latitude thermospheric 
dynamics. 

2.4. Polar Ionosphere And Polar Wind 

During the 80's it was realized empirically that the 
magnetosphere under IMF Bz northward conditions was 
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different from that under Bz southward conditions. Nowhere 
was this more clearly manifested than in the polar cap 
ionosphere. Similarly, the re-evaluation of the origin of the 
magnetospheric plasma, i.e., solar or terrestrial, has led to 
new work on the heavy ion polar wind, its transport, and its 
energization. Again, the high latitude ionosphere is a key 
region for this work. 

Our understanding of the high latitude convection under 
southward IMF conditions is complex, but well developed. 
Moses et al., [1987] considered southward IMF conditions and 
studied the plasma flow into the polar cap and its By 
dependence. Using DE-2 drift data, Moses et al., [1988] 
deduced that for about one third of the data multiple flow entry 
regions existed. Friis-Christensen et al., [1988] observed 
traveling convection vortices in the polar cleft region using a 
magnetometer network. However, these traveling vortices 
were inconsistent with present-day FTE current-closure 
models. Hairston and Heelis [1990] developed an empirical 
convection model for southward IMF conditions based upon 
DE-2 plasma drift data. Rich and Maynard [1989] developed 
analytical functions to represent the electric fields of Heppner 
and Maynard [1987]. An alternative procedure for deducing the 
electrodynamics in the high latitude ionosphere self 
consistently was developed by Richmond and Kamide [1988] 
and Richmond et al., [1990]. By combining conductivity, 
current, and electric field data, they produced the equivalent 
global current and electric field patterns which best satisfied a 
generalized Ohms law. At present this AMIE technique is 
limited by the sparsity of input data. Ahn et al., [1989] 
showed that the ionospheric conductivity dependence upon 
both the cross-polar cap potential difference and the global 
Joule heating rate were reasonably modeled by present day 
statistical models. 

Density structures in both the E and F region were the topic 
of several studies. Such density features are more readily 
observed and tracked than their electrodynamic source 
parameters. Schunk and Sojka [1987] used the TDIM to 
compare the lifetimes of F region plasma density structures and 
found that in darkness the ratio of enhanced density to the 
background density remained constant for many 10’s of hours. 
In fact, it is only by transport that these regions convect into 
regions of solar illumination or auroral precipitation and 
hence, enhanced ionization that the structures disappear. 
Sojka and Schunk [1987] similarly showed that in a dark polar 
cap in the presence of sustained multi-cell convection, regions 
of low density could be maintained because plasma transport 
would not take the flux tubes into regions of solar illumination 
or auroral precipitation and hence, enhanced ionization. 
Maynard et al., [1990] used DE-2 plasma data to search for 
those density signatures which were predicted to occur for Bz 
northward conditions. The location of the polar cap depletion 
regions was predicted to be a strong function of the IMF By 
component. Tulunay and Grebowsky [1987] compared winter 
electron density signatures in the topside ionosphere for the 
northern and southern polar regions using Ariel plasma data. 
They found systematic differences, with the southern 
hemisphere having the lowest densities. Hoegy and Benson 
[1988] compared plasma data from the ISIS and DE-2 satellites 
to study the extent and convection control of high latitude 
plasma features. Under dynamic conditions when the polar cap 
electric field was highly structured, the plasma density 
depended heavily upon this structure. Sojka and Schunk 
[1988] modeled the polar cap F region with the TDIM for these 


conditions and found that significant density depletions were 
present in the regions of strong electric Fields even if these 
enhanced fields were quite local. Heelis and Vickrey [1990] 
showed that plasma structures created in the F region produce 
plasma images in the E region through field aligned 
electrostatics. They showed how the scaling depended on the 
altitude of the source structure. 

Basu et al., [1990] used plasma and electrodynamic 
observations from DE-2 to study plasma instabilities formed at 
the edges of ionospheric patches. They contrasted the gradient 
drift with the velocity shear process, which is usually found in 
association with auroral zone density blobs. Dawn-dusk 
versus noon-midnight asymmetries were observed for the 
instabilities associated with the patches. Weber et al., [1989] 
found instability asymmetries across a polar cap sun aligned 
arc. They inferred that on the dusk side the interchange 
process was present, while on the dawn side it was the velocity 
shear instability. The effect of density irregularities was 
studied by James et al., [1990] using radio wave reception on 
DE-1 as the satellite crossed the plasma heater facility at 
Tromso. The irregularities caused both spectral broadening 
and long delays in the signals. Satyanarayana et al., [1987] 
analyzed the effects of Coulomb collisions and electron 
parallel dynamics on the transverse Kelvin-Helmholtz 
instability with the high latitude ionospheric plasma in mind. 
The parallel dynamics were found to simulate ’'buoyancy” in 
their analysis and to stabilize the instability. The effect of 
plasma turbulence in the D region associated with water cluster 
ions was analyzed by Kelly et al., [1987]. Based on 50 MHz 
radar backscatter measurements, they predicted how the 
EISCAT facility should be able to verify the water cluster ion 
theory. 

Above the polar ionosphere is a region containing 
ionospheric and magnetospheric plasma, the ionospheric 
plasma component being the polar wind. The outflow of 
plasma is controlled by various magnetospheric, ionospheric, 
and solar processes. Gombosi and Killeen [1987] considered 
how thermospheric dynamics, driven by heating, modified the 
polar wind flow. Barakat and Schunk [1987] modeled the solar 
cycle, seasonal, and magnetic activity dependence of this 
outflow. Ganguli et al., [1987] and Demars and Schunk [1989] 
developed 16 moment transport models of the polar wind. 
Gombosi and Schunk [1988] compared hydrodynamic and 
collisionless polar wind formulations for plasma expansion. 
Singh et al., [1989] studied enhanced 0+ outflows caused by 
high altitude plasma cavities, and Gombosi and Nagy [1989] 
considered the effects of return currents on the polar wind 
outflow. Cannata and Gombosi [1989] contrasted the 
composition of the polar wind between solar minimum and 
solar maximum. In their study, the solar cycle dependence was 
introduced as an underlying thermospheric dependence. This 
compositional dependence was found to be consistent with a 
vast body of satellite plasma composition data analyzed by 
Grebowsky et al., [1990]. Grebowsky et al., considered the 
altitude dependence of plasma composition in the noon +, 3 
hours sector of the polar cap for different seasonal, solar cycle, 
and magnetic activity conditions. Using DE-1 composition 
data, Chen et al., [1990] obtained the distribution of 0 + , H+, 
and He + in upward flowing beams in the polar cap. The 
horizontal transport of such beams leads to the effect called the 
cleft ion fountain. Using a three-dimensional particle code, 
Del court et al., [1989b] showed that under storm conditions, 
polar wind protons became a source for ion precipitation. 
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Horwitz [1987] showed that in describing this process, two- 
dimensional kinetic models produced different results from 
simpler one-dimensional kinetic models. Schunk and Sojka 
[1989] developed a fully three-dimensional time dependent 
model of the F region and polar wind. In this model, the polar 
wind was coupled to the ionospheric plasma and the polar wind 
was driven by temporal and spatial changes in the ionosphere. 
With such changes, transient large scale ion upflows and 
downflows occurred in the polar wind. Swift [1990] developed 
a one-dimensional fluid model to describe the ionosphere- 
polar wind from 100km to many Re. Using this model, he 
simulated changing magnetospheric flux tube boundary 
conditions, i.e., open and closed field lines. Wilson et al., 
[1990] developed a new kinetic polar plasma outflow model. 

Although ion composition measurements of the energized 
ionospheric ions as they flow into the magnetosphere are 
available, systematic observations of the thermal polar wind 
still need to be made, Sagawa et al., [1987] used DE-1 ion data 
to study the flow of 0 + and H + at L< 6 and found it to be bi- 
directional along the field lines, implying an ionospheric 
source, while Collin et al., [1988] used the same satellite data 
set to map the morphology of the He+ ionospheric outflows. 
Fuselier et al., [1989] observed two distinct 0 + populations in 
the dayside low latitude boundary layer with the AMPTE/charge 
composition Explorer. These two populations were a 
convected 0 + population from the outer magnetosphere and an 
injected 0 + population from the high latitude ionosphere. 

The question of whether the solar wind or the ionosphere is 
the principal source of magnetospheric plasma has been 
extensively investigated during the past decade. More 
recently, the idea of the ionosphere as the primary source has 
gained respect and has been quantified by a variety of studies, 
although, even today, only a few isolated observations of the 
thermal polar wind have been made. The research focuses on 
the suprathermal and energetic ion species that are clearly of 
ionospheric origin, on how they are heated and subsequently 
transported into the magnetosphere. Based upon ion 
composition measurements in the topside ionosphere- 
magnetosphere, Chappell et al., [1987] demonstrated that the 
ionosphere is a plausible source for all the plasma in the 
magnetosphere. The adiabatic motion of these ions into the 
magnetosphere was studied with a three-dimensional model up 
to a distance of 17 Re by Delcourt et al., [1989a]. In the 
absence of direct measurements of the thermal polar wind, 
identifying the location and nature of the cold plasma heating 
mechanisms is extremely difficult. Tsunoda et al., [1989] 
empirically correlated the ion fountain with the cusp while 
Collin et al., [1987] found a strong solar cycle dependence in 
the energized ion outflows. Both these results indicated that 
the initial ion heating occurred in the ionosphere. The ion 
heating mechanism is still unknown although wave particle 
mechanisms have been studied [Winglee et al., 1987; Bergman 
et al., 1988; Johnson et al., 1989; Winglee et al., 1989; Crew 
et al., 1990; Schriver et al., 1990]. 

The stability of these outflowing beams is of concern as 
they are highly non-Maxwellian and anisotropic. Barakat and 
Schunk [1987] considered the ion parallel to perpendicular 
temperature anisotropy as well as the ion to electron 
temperature ratio in the excitation of electrostatic waves. 
They found that the polar wind was very stable to such wave 
excitation and defined the temperature characteristics of this 


region of stability. The effect of convecting H+ beams on the 
stability of the polar wind was studied by Barakat and Schunk 
[1989]. Their analysis showed that the plasma was 
destabilized for relative drift energies of 1 eV and that for 
large electron temperatures the parallel propagating ion/ion 
plasma-acoustic instability was triggered first, while for lower 
values of Te the obliquely propagating ion/ion cyclotron 
instability was triggered first. 

Under northward IMF conditions the polar region becomes 
aurorally active with the presence of sun aligned arcs. 
Associated with these arcs are field aligned electrodynamic 
systems which close in the ionosphere. Several studies of 
complementary observations were carried out to deduce both 
the electrodynamics and the sun aligned arc dynamics. 
Robinson et al., [1987] combined Sondrestrom Fjord radar data 
and all sky images of a sun aligned arc to study its westward 
drift and transition into the midnight sector of the oval. 
Carlson et al., [1988] used all sky images combined with DE-2 
data and found that the arc electrodynamics were consistent 
with an electric field negative divergence model. Mende et al., 
[1988] used all sky imaging TV, the Sondrestrom Fjord radar, 
and DMSP data to follow a group of sun aligned arcs which, in 
this case, were drifting toward dawn from the midnight sector. 
They found the plasma drift inside the arc was mainly anti- 
sunward, which was also found by Robinson et al., [1987], and 
that both sunward and anti-sunward drift was found outside the 
arcs. Weber et al., [1989] launched a rocket from Sondrestrom 
Fjord which crossed a sun aligned arc and made electrodynamic 
and plasma measurements. These measurements were supported 
by aircraft all sky camera data. This arc was drifting toward 
dusk at a speed consistent with the overall electric field, 
however, inside the arc the plasma was convecting sunward. 
Inside the arc, the F region plasma showed a slight decrease 
from that outside the arc. Zanetti et al., [1990] used Viking 
magnetic and electric field data with auroral images to study the 
polar region for northward IMF conditions. They tested the 
two hypotheses for the generation of polar aurora; namely, 
converging electric fields and bifurcation of the magnetotail 
with plasma convection and found that the converging electric 
field hypothesis (V ,E < O) adequately explained the polar cap 
aurora. During the major magnetic storm on February 7, 1986, 
the IMF Bz component turned northward for several hours and 
during this time an NBZ current system developed in the 
summer polar cap while an irregular pattern of field aligned 
currents developed in the winter hemisphere [Rich et al., 
1990]. Polar cap arcs were observed but were not correlated 
with the location of the large scale field aligned currents. Chiu 
[1989] proposed that the sun aligned arcs were driven by a 
negative divergence of the electric fields, a velocity shear. In 
this model each velocity shear structure in the ionosphere 
drove a pair of field aligned currents, with the upward currents 
associated with the downward acceleration of polar rain 
electrons. 


3.1. Ionospheric Elecirodynamic Inputs 

A large fraction of magnetospheric electrodynamics links 
with the ionosphere where it can be studied in depth. A major 
part of the ionospheres ionization and energy come from the 
magnetospheric electrodynamic systems. For a complete 
analysis of the ionosphere, it is necessary to understand how 
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magnetospheric electrodynamics maps into the ionosphere; 
such a description is as yet not available. However, the last 
four years have seen significant advances in this area of 
research. Statistical magnetospheric convection models now 
include an IMF dependence, although there is no agreement on 
the convection pattern resulting from a Northward IMF 
(Heppner and Maynard 1987, Crooker 1988, Moses et al, 
1988, 1989, Lu et al., 1989). Advancements have also been 
made in describing current/convection/conductivity 
dependence upon the IMF orientation [Min-Yun et al., 1987, 
Rasmussen and Schunk 1987, 1988]. Rich and Gussenhoven 
[1987] and Hoffman et al., [1988] observed that Region 1 and 
2 current systems were virtually absent under very quiet 
geomagnetic conditions. This quiet case situation does not 
appear in present day models. Also, the statistical 
convection, current, and conductivity models do not reproduce 
the specific day-to-day variability observed in the ionosphere. 
Richmond et al., [1988] applied the AMIE procedure to infer 
the high latitude electric fields and currents during a specific 
study period. The use of global scale auroral images to infer 
the instantaneous conductivity distribution has been 
extensively developed [Marklund et al., 1987, 1988, Ahn et 
al., 1989, Rees et al., 1988]. 

Local time dependencies in the colocation of conductivity 
and convection boundaries and the subsequent effect this has 
on Joule heating was pointed out by Kamide and Richmond 
[1987], and to some extent were quantified by Foster [1987], 
while Heelis and Coley [1988] derived both global and local 
Joule heating rates based on in situ DE-2 observations. 
Coordinated observations of multiple parameters associated 
with the magnetospheric inputs to the ionosphere in the cusp- 
polar cap region were made by Basinska et al., [1987], Baker 
et al., [1989], and Niciejewski et al., [1989]. These studies 
indicate why the phenomena associated with fine structure need 
comprehensive multiparameter observation to constrain the 
physical interpretation. The latter is an example of how the 
CEDAR cooperative programs can be utilized to maximize 
constraints on the physical interpretation of ionospheric 
observations. 

Earle and Kelley [1987], using incoherent scatter radar 
data, carried out a spectral analysis of ionospheric electric 
fields at high latitudes and their penetration to the equator. 
They found a peak penetration to low latitudes for fluctuations 
in the 3-5 hour period range and noted that the zonal electric 
field components penetrated with little or no attenuation. 
Richmond and Roble [1987] used the NCAR TGCM to study the 
electrodynamic effects associated with thermospheric winds 
and reproduced the general pattern of the Sq current system. 
Aggson et al., [1987] and Ganguly et al., [1987] empirically 
studied the equatorial ionospheric convection using DE-2 data 
and Arecibo observations, respectively. The use of coherent 
HF radars to map the plasma convection in the F region was 
discussed by Ruohoniemi et al., [1989]. Optically tracking F 
region cavities generated by high power ionospheric heaters 
to deduce plasma convection was demonstrated by Bernhardt et 
al., [1989]. 

Extensive conjugate observations using satellites in both 
polar regions have quantified the extent to which polar cap 
arcs and the cusp are conjugate [Mizera et al., 1987; Candidi 
and Meng 1988; Carbary and Meng 1988]. Newell and Meng 
[1988] reported the observation of an ion polar rain event 
which was significantly larger than any previously reported, 
while Sergeev et al., [1990] found transient plasma injection 


events which were not linked with substorms. 

The role played by the ionosphere in the transmission of 
VLF signals, whistlers, and lightning continues to be a topic 
where progress is being made. Inan et al., [1989] 
quantitatively verified the predictions of lightning-induced 
electron precipitation (LEP) by observing the gyroresonant 
whistler particle interaction on an L=2.24 field line along with 
a burst of precipitation. A longitudinal dependence of electron 
precipitation at L=4 associated with the south Atlantic 
anomaly was observed by Bering et al., [1988]. Neubert et 
al., [1987] proposed the resonance of coherent whistler mode 
waves with energetic electrons in the topside ionosphere as a 
source of Fast Trimpi events. The modification of sub- 
ionospheric VLF signals by LEP ionization changes in the D 
region was observed and statistically modeled [Carpenter and 
Inan 1987, Inan and Carpenter 1987, Inan et al., 1988]. 
Baginski et al., [1988] used a model of the lightning return 
stroke to show that it can be a major component of the 
transient field and energy coupled to the ionosphere. The role 
of mid-latitude precipitation in modifying the geoelectric field 
was studied by Sheldon et al., [1988] who found that the 
precipitation produces a diurnal modulation of the geoelectric 
field. 

3.2. Plasma Waves And Irregularities 

The ionosphere, apart from having a diverse collection of 
particle species, harbors a zoo of naturally present wave- 
particle processes. Many of the research topics already 
reviewed have covered some of these processes. However, 
other studies have been conducted which are of a more general 
plasma physics nature, or which do not readily fit within the 
constraints of the previous sections. These are reviewed in 
this section. Much of this work involves solving non-linear 
systems of equations whose boundary conditions are not well 
understood; nonetheless, significant progress has been made 
on several of these wave-particle problems. 

Maggs [1989] considered the nonlinear spatial evolution 
of an auroral electron beam and of beam generated electrostatic 
whistler noise as the beam propagates into the atmosphere and 
found that the wave intensities were sufficient to modify the 
ionospheric density profile. Seyler [1988] considered the 
three-dimensional nonlinear evolution of kinetic Alfven 
waves due to shear flow and collisionless tearing instability 
under boundary conditions similar to those of a discrete auroral 
breakup arc. Mode conversion to fast magnetosonic waves by 
anisotropic ionospheric currents from shear incident Alfven 
waves was modeled by Fujita and Tamao [1988] and Fujita 
[1988]. They considered the dispersion characteristics, loss 
mechanisms, and the ground magnetic signature of this mode 
conversion mechanism. The heating of ions in a 
multicomponent plasma was found to occur even when the 
electron drift speed was too small to excite the lower hybrid 
instability. In this case, Ashour-Abdalla et al., [1987] found 
that an ion-ion hybrid mode became unstable and heating of 
both the light (bulk heating) and heavy (tail heating) ions 
occurred [Schriver and Ashour-Abdalla 1988]. Suscynsky et 
al, [1989] carried out a laboratory two-ion plasma experiment 
to generate EIC waves. They used K + for the light ion and Cs + 
for the heavy ion and found that the EIC wave frequencies 
depended on the density ratio of the ions. Benson and Wong 
[1987] used ISIS-1 observations to show that the cyclotron 
maser instability generates most of the auroral radio 
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emissions, and they predicted the characteristics of these radio 
emissions for ground observations. A ground-based search for 
these auroral radio emissions found them to occur within ± 30 
minutes of magnetic midnight but that they were poorly 
correlated with geomagnetic activity [Benson et al, 1988a]. 
Using DE-1 data, Benson et al, [1988b] found ordinary mode 
AKR fine structure; however, the present day theories used to 
describe the intense right-hand extraordinary mode AKR do not 
generate this left-hand ordinary mode/AKR. Keskinen et al., 
[1988] considered the non-linear evolution of the Kelvin- 
Helmholtz instability caused by velocity-sheared plasma 
flows. Specifically, they studied the role of ion-neutral 
collisions and compared their results with DE observations. 

Plasma wave interactions were the topic of several studies; 
De Groot and Mizuno [1989] considered the effect of the ion 
acoustic decay instability on microwave-plasma interactions. 
Yu et al., [1987] showed that in an ionospheric plasma a large 
amplitude electromagnetic wave can parametrically excite low- 
frequency electrostatic modified electron acoustic waves; 
Vampola [1987]; Gail et al., [1990]; Gail and Inan [1990]; 
Datlow and Imhof [1990]; Rosenberg et al., [1990] and others 
identified in Section 2.3 considered the role of VLF waves in 
the precipitation of ring current particles; Groves et al., 
[1988] observed intense hiss from an auroral rocket flight and 
compared their observations favorably with the predictions of 
a convective beam amplification hiss model. A variety of 
plasma wave phenomena can be generated using the HF 
sounders and heating facilities; Osherovich [1987] discussed 
the physics of the diffuse plasma resonances as stimulated by 
HF sounders; the resonant absorption of HF waves by the 
electrons to produce a "bump in the tail" distribution was 
studied by Shoucri et al., [1987] and Villalon [1989]; while 
Hansen et al., [1989] refined the calculations of how deep the 
associated electron density depletions (cavities) can become. 
Swartz et_ al., [1988] discovered that the interpretation of 
topside incoherent scatter echoes can be heavily aliased if 
velocity shears are present in the plasma. This effect tends to 
reduce the apparent electron temperature while enhancing the 
ion temperature, even to the extent of wrongly deducing that 
Te < Ti, Lockwood et al., [1988] and Provjdakes et al., [1988] 
considered the role of plasma temperature in the F and E 
regions, respectively, in wave plasma interaction using 
ground-based sounding techniques. 

Basu et al., [1987] used an HF heater facility to study the 
growth and decay rates of ionospheric scintillations associated 
with 750m scale irregularities. Villain et al., [1987] used 
another HF radar to study E region decametric irregularities 
produced by oblique electron streaming. The HF heater facility 
was also used to produce Field aligned irregularities which had 
spread F echoes associated with them; in this case, a thermal 
instability generated the density irregularities [Kuo and Djuth, 
1988]. Bell and Ngo [1990] showed that electrostatic lower 
hybrid waves are excited by whistler mode waves scattering 
from field aligned plasma density irregularities. A generalized 
Rayleigh-Taylor instability model was used to study the 
collisional interchange instability associated with equatorial 
irregularities and the results were successfully compared with 
observations [Zargham and Seyler 1987; Kelley et al., 1987], 
Zargham and Seyler [1989] included ion inertia and ion-neutral 
collisions in a two-dimensional fluid model of ionospheric 
interchange instability. Using this model, they computed the 
density structure power-spectra for several forms of plasma 
irregularities: rising plasma bubbles, unstable plasma strata, 


and driven plasma irregularities. Lee et al., [1989] modeled the 
effect of both rod-like and sheet-like ionospheric irregularities 
on Faraday polarization fluctuations. They found that the 
sheet-like irregularities had the most severe effect on radio 
waves; such irregularities being found in both the equatorial 
anomaly and auroral regions. The fading of HF waves returned 
from the F region when the receiver is close to the transmitter 
was studied by Booker and Tao [1987]. Franke and Liu [1987] 
used fading of saturated VHF amplitude scintillations to deduce 
the characteristics of equatorial plasma bubbles. 

Lakhina et al., [1990] considered the stability criterion of 
the ballooning mode instability at the plasmapause; Hollweg 
[1990] studied the resonance absorption of propagating fast 
waves in a plasma density gradient; while Sentman [1989, 
1990] studied the polarization and mode splitting in Schumann 
resonance excitation as well as the resonance phenomena in a 
two-scale-height ionosphere. Liao et al., [1989] proposed the 
parametric instability excited by lightning-induced whistler 
waves as a cause of explosive spread F. 

3.3. Active Experiments 

The ionosphere is a readily accessible "unbounded" plasma 
suitable for basic plasma physics research. Some of this 
research also highlights certain ionospheric physics 
problems. This section focuses upon this subset of research. 

Two major high power HF heater facilities at Arecibo, 
Puerto Rico and at Tromso, Norway have been used extensively 
to modify the upper atmosphere. The heating of ionospheric 
plasma, specifically the electrons, was reported by Fejer and 
Sulzer [1987], Bernhardt and Duncan [1987], Djuth et al., 
[1987b] and Inan [1990], while Shoucri et al., [1987] 
described how the electron distribution function was modified 
by the HF waves through a resonant electrostatic field 
interaction. Newman et al., [1988] pointed out that under 
certain solar/seasonal conditions the apparent anomalous 
heating was due to very low electron cooling rates. 
Osherovich [1987, 1989] developed theoretical descriptions 
of the stimulated diffuse plasma resonances observed during HF 
heating experiments. Chaturvedi et al., [1987] described the 
production and control of ion-cyclotron instabilities at high 
latitudes using high power radio waves. Hansen et al., [1989] 
studied theoretically the lifetime of density depletions created 
by the heating experiments, while Bernhardt et al, [1989a] 
showed how these density depletions could be used as a tracer 
for ionospheric convection and also how the associated 
optical signatures could be observed and tracked [Bernhard et 
al, 1989b], Associated with these density gradients are waves 
and irregularities; Bernhardt [1988] showed how the E x B 
gradient drift instability produced irregularities on the upwind 
edge of a density depletion. Using the HF heater to create 
density depletions, Basu et al, [1987] was able to observe the 
e-folding growth and decay rates of F region irregularities of 
the 750m scale size while Noble et al, [1987] made a similar 
study of the 1 to 7m E region irregularities. The mechanism by 
which the HF experiments are able to stimulate enhanced 
plasma line emission and the "overshoot" aspects are still 
unknown. Observations of the enhanced plasma line and the 
"overshoot" were made by Djuth et al, [1987a] and Djuth and 
Gonzales [1988], and the role of Langmuir waves in the 
overshoot process was studied theoretically by Kuo et al, 
[1987]; however, Sulzer et al, [1989] showed that significant 
aspects of the phenomena still remain unexplained. 
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The artificial generation of F region bubble irregularities 
was reported by Klobuchar and Abdu (1989). Rockets were 
launched into the bottomside F layer from Natal, Brazil as part 
of the BIME program. These payloads released FfeO and CO2 to 
create plasma depletions. The subsequent motion of the 
depletions and evolution of plasma irregularities provided a 
good test of equatorial bubble irregularity theories. 

In-situ generation of whistler mode waves by an electron 
beam generator and a remote wave detector on board the space 
shuttle was the focus of both observational and theoretical 
reports. On both the space shuttle STS -3 and SpaceLab-2 
flights the electron generator (FPEG) and a plasma monitor 
(PDP), respectively produced and observed enhanced VLF 
emissions [Farrell et al., 1988; Frank et al., 1989; Reeves et 
al., 1988, 1990]. Goerke et al., [1990] observed VLF 
emissions generated by an electron beam gun on a Nike Black 
Brant V sounding rocket; similarly, Winglee and Kellogg 
[1990] reported electromagnetic wave emissions generated by 
beams, and Olsen et al., [1990] studied plasma waves generated 
by ion beams on the SCATHA satellite. The role of coherent 
Cerenkov radiation driven by bunching of the electron beam 
was developed theoretically by Farrell et al., [1989] and 
Pickett et al., [1989]. Nishikawa et al., [1989] modeled the 
whistler mode excitation using a three-dimensional partially 
magnetostatic code and subsequently identified the beam 
instability responsible for the excitation. In a theoretical 
study, Okuda and Ashour-Abdalla [1988] showed that 
suprathermal electron beams could also drive ion-acoustic 
waves in the ionosphere which, in turn, leads to ion heating. 
Water dumps from the space shuttle that produce a host of 
plasma irregularities in the few Hz to KHz range were studied 
by Pickett et al., [1989], who considered both ion-plasma and 
Ott-Farley instabilities as sources of plasma density 
irregularities. 

In addition to these active experiments in ionospheric 
physics, several other categories of ionospheric modification 
such as beams in the ionosphere, chemicals in the ionosphere, 
and ionosphere-vehicle interactions were reported. As these 
studies apply primarily to plasma physics, only a 
categorization and references are given below. Various aspects 
of the plasma vehicle interaction, i.e., wake, charging, 
sheath, were studied (Lai et al., [1987]; Samir et al., [1987]; 
Okuda and Berchem [1988]; Singh and Kwang [1988]; 
Hastings and Blandino [1989]; Katz et al.,[1989]; Kaufman et 
al., [1989]; Ma and Schunk [1989]; Murphy and Katz [1989]; 
Murphy et al., [1989]; Myers et al., [1989]; Samir et al., 
[1989]; Senbetu and Henley [1989]; Tribble et al., [1989] 
Kellogg et al., [1990]; Neubert et al., [1990a, b]; Svenses et 
al., [1990]). Chemical releases leading to plasma 
modifications and neutral clouds were reported by Bernhardt 
[1987]; Bernhardt et al., [1987]; Grebowsky et al., [1987a, b]; 
Hunton et al., [1987]; Lai et al., [1988]; Eccles et al., [1989]; 
Lai and Murad [1989]; Paterson and Frank [1989]; Elgin et al., 
[1990]; Gilchrist et al., [1990]; Jacobson et al., [1987]; Ma 
and Schunk [1990]; Pike et al., [1990]. The physics of beams 
in ionospheric plasmas was extensively reported upon (Cai et 
al., [1987]; Chung [1987]; Harker and Banks [1987]; Parish et 
al., [1987]; Pritchett and Winglee [1987]; Banks and Raitt 
[1988]; Hudson and Roth [1988]; Pollock et al., [1988]; 
Reeves et al., [1988]; Abe et al., [1988]; Erlandson et al., 
[1989]; Frank et al., [1989]; Heppner et al., [1989]; Neubert 
and Harker [1988]; Stenzel and Urrutia [1989]; Hallinan et al., 
[1990]; Mandell et al., [1990]; Pritchett [1990]; Scales and 


Kintner [1990 a,b]; Winglee [1990]). Additional HF plasma 
heating results were described by Isham et al., [1987]; Fejer et 
al, [1989]; Kuo and Lee [1989]; Thide et al., [1989]; Kuo and 
Lee [1990]; Milikh [1990]. The critical velocity ionization 
phenomenon was studied by Lai et al., [1990 a,b]; McNeil et 
al., [1990]; Person et al., [1990]; Stenbaek-Nielsen et al., 
[1990] and others already listed above. 

3.4. Ionospheric Forecasting 

Since ionospheric physics is regarded by many as a mature 
research field, one measure of its maturity is the quality of the 
community’s modeling and forecasting ability. The primary 
parameters to be modeled and predicted are the plasma density, 
its temperature and composition, and the waves generated in 
the presence of anisotropic plasma distributions. Models 
exist for each of the regions described in Section 2; however, 
most of these models are limited, either inherently by physical 
assumptions or by inadequately defined boundary conditions. 
In this quadrennium, improvements in modeling have occurred 
in several areas. These are discussed in the following 
paragraphs. 

The RIMS instrument on the DE-1 satellite produced 
measurements that led to new insights into the plasmasphere s 
exotic, but normal status. The role of minor ions, their solar 
cycle, and storm dependencies have been incorporated into 
recent models of the plasmasphere [Chandler et al., 1987, 
Newberry et al„ 1989]. The question of where refilling begins 
along a plasmaspheric flux tube was addressed by Rasmussen 
and Schunk [1988] using a new multistream hydrodynamic 
model. They confirmed earlier findings that it is at the equator 
that refilling begins, but showed that in the initial refilling 
period (4 hours) the flow characteristics were substantially 
different from those previously reported. Another example 
showing the limitations of present day models is associated 
with electron temperature modeling and global scale coupling. 
Brace et al., [1988] used both DE-l and DE-2 plasma 
measurements to correlate the location of plasmaspheric and 
ionospheric enhanced electron temperatures and found that in 
the evening bulge region they were not correlated. This bulge 
region is very dependent upon magnetic activity, and, 
therefore, this discrepancy is readily explained in a qualitative 
way. Present day models do not reproduce the dynamics of 
such ionospheric plasmaspheric temperature structures. 

Even greater difficulty is encountered in describing and 
modeling observations at higher latitudes. Persoon et al., 
[1988] found that in the topside ionosphere (polar and auroral 
flux tubes) marked plasma depletions and structuring exist. In 
such regions, the cold plasma is susceptible to acceleration by 
a wide range of auroral processes which are not included in 
present day models, e.g., Demars and Schunk [1987], Schunk 
and Sojka [1987, 1989]. Even when auroral systems are 
included, i.e., return currents as by Schunk et al., [1987], the 
work is still of a parametric nature. It is difficult to include 
micro scale phenomena in large scale models. St. Maurice 
[1987] outlined how, in the case of E region electron heating 
in plasma irregularities, the micro scale heating rate could be 
parameterized on a larger scale. Rasmussen et al., [1988] and 
Richards et al., [1989] showed the difficulty of using present 
day models and extensive observations to describe case studies 
in the sub- auroral/mid-latitude F region. In both case studies, 
significant, unexplained differences remain, although ex- 
tensive validation of the model inputs was carried out. 
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Buonsanto [1989] compared Millstone Hill incoherent scatter 
radar data with the International Reference Ionosphere (IRI) 
model for different solar, seasonal, and diurnal conditions. 
The topside densities were significantly larger in the model 
while the observed Te variability was not present in the IRI 
model. The observed Ti diurnal variation was present in the 
model but with a small amplitude. Sica et al M [1990] used 
extensive data sets from many locations to test the TDIM mid 
-latitude model. A similar test was carried out by Sojka et a!., 
[1988] using a single station data base over an entire solar 
cycle. In both these cases, since the inputs were not 
adequately constrained, agreement could be obtained in several 
ways. Often the inputs are still inadequately understood as is 
the case for the cleft plasma heating source. Fontheim et al., 
[1987] showed that elevated electron temperatures (up to 
10,000° K) were independent of solar EUV, but were coupled to 
precipitation and, hence, concluded that particle heating was 
the source mechanism. However, to include such a mechanism 
in present day ionospheric models requires much improved 
knowledge or parameterization of magnetospheric cleft 
precipitation. 

On a global scale only climatological ionospheric 
modeling is possible at present. Several studies were 
conducted and models developed to improve our global scale 
climatological knowledge. Roble et al., [1987] carried out a 
global mean coupled ionosphere-thermosphere study; 
Anderson et al., [1989] developed a fully analytic mid- and low 
- latitude model, while Sojka and Schunk [1989] contrasted F 
region seasonal global effects. In each case, although new 
capability was achieved, forecasting capability was only at a 
climatological level. To extend our capability beyond this 
climatological status, dynamics associated with the regions 
coupled to the ionosphere, namely, the magnetosphere and 
thermosphere, are needed. Knipp et al., [1989] used extensive 
data sets and the AMIE technique to produce storm dependent 
electrodynamic patterns for a specific study period. Sojka et 
al., [1989] used DE-1 auroral images to define auroral dynamics 
for an ionospheric simulation. Since the associated 
electrodynamics were not available, a "parametric" study 
resulted rather than an actual ionospheric simulation. 

The accuracy of our ionospheric physics models is still far 
from being predictive. At best, climatological predictive 
modeling is possible for certain phenomena, however, these 
are mostly of a semi-empirical nature [Schunk and Sojka, 
1988]. At the present time, complete understanding of the 
phases and morphology of an ionospheric storm illude us. 
Part of the reason for this fs the lack of understanding of how 
inputs to the ionosphere behave (magnetospheric and 
thermospheric dynamics), the difficulty of including small 
-scale size variability in global scale simulations, and of 
understanding how certain processes operate (auroral plasma 
heating, irregularities, return currents, etc.). All of this can be 
summarized by the need to include ionospheric weather rather 
than only deal with climatological trends. During this 
quadrennial period, preliminary studies have examined the 
spatial and temporal e ffects of ionospheric weather on 
climatological models. These studies have also highlighted 
the inability to tweak climatological model inputs (statistical 
models for auroral precipitation, plasma convection and 
neutral atmosphere) to give weather. Rich and Maynard [1989] 
carried out a critique of statistical representations of the 
electrodynamic input to the ionosphere and raised some of 
these issues. 


3.5. COUPUNG THE IONOSPHERE WITH 
Other Regions 

Both theoretically and experimentally the community 
continues to make progress in coupling other regions to the 
ionosphere. The characteristics of hydromagnetic waves in 
the magnetospheric cavity are, to a large extent, determined by 
the ionosphere. Crowley et al., [1987] observed Pc-5 
pulsations using ground based magnetometers and EISCAT and 
were able to deduce that the Pc-5 damping was controlled by 
the damping of the cavity mode rather than the energy 
dissipation from the resonance field line. Lysak [1988] 
showed theoretically how bursts of irregular pulsations (PiB) 
were associated with a resonant cavity in the topside 
ionosphere where the density decreased exponentially with 
height, and Kivelson and Southwood [1988] reconsidered the 
role of the ionosphere in the theory of hydromagnetic surface 
waves in the magnetosphere. The ionosphere is also a key 
element in auroral dynamics and arc formation. Lotko et al., 
[1987] introduced the ionosphere into a two-dimensional 
model of perpendicular flow dynamics in boundary layers in 
the magnetospheric equatorial plane coupled with the 
ionosphere and found that the ionosphere controlled the 
disturbances at different scale sizes. Kan et al., [1988] found 
that the co-location of the magnetospheric convection 
reversal with a conductance gradient in diffuse aurora 
determined the stability and onset of substorms. Zhu and Kan 
[1987] showed that the dynamics of westward traveling surges 
in a coupled magnetospheric-ionospheric system was 
controlled by the recombination rate in the auroral 
ionosphere. Mapping electric Fields from the ionosphere into 
the magnetosphere during disturbed times in the vicinity of 
westward traveling surges was studied by Kaufmann and Larson 
[1989]. Chiu [1989] argued that the magnetospheric 
convection velocity shears in the ionosphere drove and 
supported polar cap auroral arcs without the need to bifurcate 
the magnetospheric topology with adjacent open and closed 
field lines. Coupling the polar wind from the magnetosphere 
to the polar ionosphere was accomplished theoretically by 
Schunk and Sojka [1989] with a three-dimensional time 
dependent model. This model enabled ionospheric topside 
boundary conditions to be replaced with more realistic 
physical processes in the high altitude polar magnetosphere. 

Solar x-ray fluxes can lead to photoelectrons in the 
ionosphere at energies higher than those associated with the 
solar EUV spectrum and, hence, provide the magnetosphere 
with an additional electron source [Winningham et al., 1989]. 
Transient lightning pulses transmit energy into the 
ionosphere and magnetosphere; Baginski et al., [1988] self 
consistently coupled these transient pulses with the 
ionosphere and showed that significantly more energy than 
previously expected was directed to the ionosphere. Inan and 
Carpenter [1987] showed how lightning induced electron 
precipitation caused modification of the D region conductivity 
which, in turn, modified the ionosphere’s role in VLF 
transmission. The modification of the D region by relativistic 
electrons associated with high speed solar wind streams was 
inferred by Baker et al., [1987]. 

Coupling of the ionosphere and thermosphere has 
progressed significantly. Roble et al., [1988] has developed a 
self-consistent thermosphere-ionosphere model (TIGCM) 
which combines the physical processes of each. Work 
continues in deducing self consistent ionospheric and 
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thermospheric parameters from the analysis of incoherent 
scatter radar data [Burnside et al., 1988]. Johnson and 
Luhmann [1988] were able to infer from the lack of 
geomagnetic control of upper mesopheric neutral winds that 
magnetospheric ion drag coupling was only significant above 
90km at high latitudes. Propagation of magnetospheric storm 
effects to low latitudes in both the ionosphere and 
thermosphere was the focus of an ionosonde study by Forbes et 
al., [1988]. The recovery phase of another storm was studied 
by Burns et al., [1989] to investigate the thermospheric 
composition variations. The effect of neutral winds on 
magnetospheric-ionospheric electrodynamics was studied by 
Forbes and Harel [1989]. The neutral wind control of 
electrodynamic coupling in the equatorial E and F regions was 
studied by Biondi et al., [1988]. Samson et al., [1989] used 
the Goose Bay coherent scatter radar to show that dayside 
atmospheric gravity waves originated at an ionospheric 
convection shear reversal region in the 100 to 150km altitude 
range. 

Conclusion 

Over the past four years advancements have been made in 
all areas of ionospheric research. As already pointed out, the 
ionospheric community is a very balanced research community 
and, as such, contributes to our knowledge over a diverse range 
of sub-topics. Some of the work leads to greater confusion, 
again due to the complexity of the phenomena and paucity of 
observation. In several areas, significant new results were 
published during this reporting period. Some highlights are 
listed below. 

• Global images of auroral emissions have been used to 
infer the temporal and spatial pattern of auroral conductances 
and electron precipitation. However, the single pixel 
representation of the fine-scale auroral structure is still in 
debate. This is yet another example of how phenomena on 
different scale sizes need to be interrelated. 

• In the same vein, the need for instantaneous global 
"images" of the magnetospheric convection pattern has been 
realized and, although it is in its infancy, the AMIE technique 
has demonstrated its usefulness. The community has a great 
need for global scale time-dependent observation of the 
convection pattern. 

• Storm dynamics in the ionosphere is a documented, but 
poorly understood phenomenon. However, progress has been 
made in quantifying the storm propagation from high to low 
(equatorial) latitudes. This is especially the case for the 
electric field penetration, F-layer dynamics, and coupled 
ionospheric-thermospheric effects. We still do not have an 

adequate representation of the storm input dynamics or its 

energy input profile at high latitudes. 

• Much work was focused on understanding the high 
latitude ionospheric-electrodynamics associated with 
northward and southward IMF conditions. Current systems 
(NBZ), FTE signatures, polar cap arcs, and ionospheric 
density/ temperature signatures were studied. However, the 
major question of the northward IMF convection pattern 
topology, namely "distorted 2-cell" or "4-ceH", remains 
unanswered. 

• Many plasma instability mechanisms which lead to 
ionospheric density irregularities were quantified. A better 
understanding of the initial ionospheric condition associated 
with these instabilities was obtained. The task of how to 


model the competing instability mechanisms is now in hand. 

• The interaction of small-scale plasma processes in a 
global ionospheric electrodynamic system was studied. 
Extensive complementary data sets were acquired arid 
rudimentary simulations were carried out including coupling 
into a magnetospheric environment. A need exists for the 
parameterization of these small-scale processes so that their 
impact on global scale models can be understood. 

• Active experiments were used to study growth and decay 
rates of plasma bubbles, density cavities, and irregularities 
both in the E and F regions. Such active experiments give 
quantifiable insight into an otherwise complex, if not 
apparently turbulent, phenomena. 

■ Ground or lightning generated VLF waves ducted through 
a "magnetospheric-ionospheric” system interacting with 
energetic magnetospheric ring current particles were studied 
via their ensuing precipitation into the atmosphere. 
Quantitative modeling of this coupled wave-particle system 
was successfully carried out as was the ensuing modification of 
the ducting channel by the precipitation. 

• Ionospheric-magnetospheric coupling via the transport 
of plasma into the magnetosphere was quantified both at high 
latitudes (polar wind, ion conics, and ion beams) and at mid- 
latitudes (plasmaspheric ion composition). Indeed, with the 
exception of the cold polar wind, this coupling has now been 
determined empirically although the energization of ions in 
conics and beams requires further work. 

• Coupling of the ionosphere to the thermosphere has seen 
the evolution of global scale coupled models. These models 
treat the thermospheric and ionospheric problems self 
consistently, and, hence, have the potential of better 
simulating storm dynamics, a topic of major interest to both 
communities. 

Glossary 

AE: NASA atmospheric satellite program. 

AKR : Auroral Kilometric Radiation. 

AMIE. Assimilative Mapping of Ionospheric Electro- 
dynamics, a computer package to best fit ionospheric currents, 
conductivities, and electric fields as global patterns. 

AMPTE: International solar wind-magnetosphere satellite 
program. 

Blob : Localized auroral region of enhanced plasma density in 
the F-region. 

CEDAR: Coupling Energy and Dynamics of Atmospheric 

Regions, 

CONDOR: Joint USA-Brazilian Rocket Program. 

CUPRI : Cornell University Portable Radar Interferometer. 
DE: NASA thermosphere-ionosphere-magnetosphere satel- 

lite program. 

DMSP : U.S. Air Force Defense Meterological Satellite 

Program of satellites. 

EIC: Electrostastic Ion Cyclotron plasma waves. 

EISC AT: European Incoherent SCATter radar facility. 

ETS: Equinox Transition Study of the thermosphere and 

ionosphere. 

EUV: Extreme Ultra Violet region of the solar spectrum. 
FLIP : Field-Line Interhemispheric Model of the University 
of Alabama, Huntsville. 

FTE: Flux Transfer Event. 

Fof2 : The plasma frequency associated with the F-region 
peak electron density. 
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FPEG: Fast Pulse Electron Generator flown on the Space 
Shuttle. 

GISMOS: Global Ionospheric Simultaneous Measurements 

Of Substorms. 

HILAT: USAF ionosphere-magnetosphere satellite. 

hmf2 : The altitude of the peak electron density in the F- 

region. 

IMF : Interplanetary Magnetic Field of solar origin. 

IRI : International Reference Ionosphere model. 

ISIS: Canadian ionospheric satellite program. 

LEP : Lightning-induced Electron Precipitation from the ring 
currents. 

NBZ: Northward orientated interplanetary magnetic field. 
NCAR: National Center for Atmospheric Research. 

Ntnf2\ The maximum electron density in the F-region profile. 
NOAA: National Oceanic and Atmospheric Administration. 
MITHRAS : Magnetosphere Ionosphere THermosphere 

RAdar Studies. 

MSIS: Mass Spectrometer/Incoherent Scatter atmospheric 

model. 

MU: Middle and Upper atmosphere VHF radar, Kyoto, Japan. 
Patches : Localized polar cap regions of enhanced plasma 

density in the F-region. 


PDP: Plasma Diagnostic Package flown on the Space 

Shuttle. 

RIMS: Retarding Ion Mass Spectrometer instrument on the 
DE-1 satellite. 

SAR: Stable Auroral Red arc. 

STS: NASA Space Transport System (shuttle). 

SUNDIAL: Coordinated study of the ionosphere/ 

magnetosphere. 

TDIM : Utah State University Time Dependent Ionospheric 
Model. 

TEC: Total Electron Content along a ray path through the 
ionosphere. 

TGCM: NCAR Thermospheric Global Circulation Model. 
TID: Travelling Ionospheric Disturbance. 

TIROS: NOAA weather satellite program. 

UV: Ultra Violet region of the solar spectrum. 

VHF : Very High Frequency radio wave. 

Viking: Swedish ionospheric-magnetosphere satellite. 

VLF: Very low frequency radio waves. 
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ABSTRACT 

Our understanding of the global scale magnetospheric electrodynamics system will, even after 
ISTP missions, still be only qualitative. Indeed, it is only the auroral imagers which enable us to 
have a 2-D time-dependent glimpse of the complexity of the magnetosphere. The imagers will, 
however, not be able to map the key electric Fields and currents responsible for the complex 
dynamics of the magnetosphere. A 2-D time-dependent look at the electric fields and currents is 
desperately needed. 

Alternative space measurement technologies must be pursued to achieve mapping of these key 
electrical fields. A system consisting of a constellation of 50 small satellites is described to achieve 
the "imaging" of these parameters. Each small satellite will support a single, simple electric field 
(and/or current) measurement instrument, and will use off-the-shelf technology to support the 
experiment Data sampling rates, onboard storage, and telemetry will be minimized to achieve the 
single overall scientific goal; namely, the global 2-D mapping of the electric field (and/or current) 
from altitudes of 500 to 800 km, every fifteen minutes. 

INTRODUCTION 

A major aspect of scientific research in the Thermospheric-Ionospheric-Magnetospheric (TIM) 
system has been via community wide collaborative campaigns and programs (GISMOS, CEDAR, 
WITS, CD AW, etc.); the object being to better understand the larger scale phenomena. In most 
cases die global scale is not observed, although the advent of global auroral imaging is beginning 
to change this. For most of these campaigns and programs an understanding of the global 
transport of energy, or mass, or electrodynamics is the major goal. At present, die observations 
are limited to local observations or ID satellite coverage. 

The future ISTP program will be similarly restricted. Indeed, the magnetospheric community 
identifies these 3-D problems as major questions not satisfactorily addressed by the ISTP program 
[Banks et al., 1988]. These same problems are also critical to the future GEM program [Roederer 
et al., 1988]. Both of these reports are based on community wide input 

A global distribution of simultaneous electric field observations is needed. Both hemispheres need 
to be covered to identify IMF and other hemispheric asymmetries. No global imagery of the 
electric field is currendy possible, and ground based global coverage on this scale is impractical. 
Hence, in situ measurements need to be made. A single satellite produces a ID cut and therefore 
to get global coverage a constellation i.e. a fleet of satellites is needed. In this paper we consider a 
constellation of 50 small satellites. By varying the individual orbit parameters, i.e., ascending 
node, inclination, UT of ascending node, etc., an optimized coverage is obtained. Circular orbits 
whose altitudes range from 500 to 1000 km would be sufficient Each satellite could produce as 
little as one electric field measurement every 0.1 second (i.e. 1 km. along the satellite track). 

Higher frequency data collection would not produce significantly better 2-D coverage or resolution. 




PRESENT DAY MAGNETOSPHERIC ELECTRIC FIELD STATUS 


A wide variety of empirical and semi-empirical models of the magnetospheric electric field as 
mapped into the high latitude ionosphere currently exists [Volland, 1978; Heppner, 1977; Heelis et 
al., 1982; Foster, 1983; Heppner and Maynard, 1987], However, these models are based upon 
either extensive satellite data bases or ground based incoherent scatter radar data sets. In either 
case the instantaneous or simultaneous data coverage is extremely limited, to either a satellite track 
or locally high resolution ground radar data. By binning (or sorting) the data according to 
magnetic or solar indices, data sets of improved spatial coverage are achieved. This is the basis for 
the empirical models. 

The resulting models no longer have temporal (substorm or even storm) dependencies. Even their 
gross morphologies are contested. Our community is faced with the quandary of deciding between 
"distorted two-cell" versus "four or multiple cell" descriptions of the convection electric field under 
IMF Bz northward conditions. Since both interpretations are based upon the same satellite data 
sources, namely, the Dynamics Explorers, there is great difficulty in finding "ground truth" data 
against which the models can be compared. 

In specific study periods where a global picture of the convection is needed, several methods are 
used to construct this convection pattern. One method is to scale an empirical electric field for the 
required magnetic and IMF conditions. A parameter such as the cross polar cap potential [Reiff 
and Luhmann, 1986] is computed and used to define the scaling factor. A more eleborate scheme 
is to define the polar cap potential along its boundary and then use this to define the global potential 
[Lu et al., 1989]. In either case, dynamics is introduced via the change in magnetic or solar indices 
and hence different choices of convection patterns. Figure 1 shows an example of how difficult it 
is to define this boundary potential from limited satellite coverage [from Lu et al., 1989]. Two 
orbits of DE-2 are compared with the corresponding global convection patterns. The electric fields 
in the cusp and Harang regions are very poorly defined. 

A second method is to use all the electrodynamic data available, currents, conductivities, and 
electric field data to self-consistently deduce the best global picture of the currents and electric 
fields. This technique has been pioneered by Richmond [Richmond et al., 1988] and is referred to 
as the Assimilative Mapping of Ionosphere Electrodynamics (AMIE). At present, even this 
technique is unable to produce significant results because of the lack of data. The process relies 
heavily upon a statistical electric field model when data is sparse, and consequently, it loses its self 
consistent attributes. Figure 2 shows an example of how sparse the data distribution typically is 
[from Richmond et al., 1988]. Using ground based radars and magnetometers as well as satellite 
data, the coverage is still too sparse to be meaningfully inverted into a global scale convection 
pattern. 

Although ground based coverage could be somewhat enhanced the cost of new facilities alone 
would be prohibitive, and these facilities would be unable to operate continously. At present, the 
incoherent scatter type radars give the most extensive and reliable data, but these radars can only 
operate for a few days per month. Installing one or two more such facilities would clearly not 
result in a major difference in coverage. They would produce additional, but still only local 
measurements. 


THE ELECTRIC FIELD (ELF) SATELLITE 

The design and development of these small ELF satellites has to take on novel measures to ensure 
that 50 can be procured (afforded), deployed (globally), and yet give good scientific quality data. 
This was the goal of Dr Frank Redd's* NASA sponsored Space Design Course for 1989-1990. A 
group of 18 students participated in the design of the various satellite subsytem/ mission 



components of such a 50 satellite mission. A final project report has been prepared by Redd and 
Olsen [1990]. The salient features of this study are summarized in the following subsections and 
Section 4. 

3.1 Design Considerations. The primary performance driver for this mission is the need to measure 
the attitude of each satellite very accurately while making the 3 axis electric field measurements, hi 
addition, it is necessary to know the induced electric field generated by the satellite as it crosses the 
magnetic field lines. This must then be factored out of the measurements. The electric field booms 
must rotate in order to balance photoelectric effects and aid in the measure of the V X B bias. 

3.2 Electric Field Sensing System. The 3 axis system will consist of 3 orthogonal sets of insulated 
booms with conductive spheres attached to the ends as shown in Figure 3. The electrical potential 
across each boom will be measured. Photoelectric charge build up on the booms is minimized by 
having the satellite rotate. To achieve the desired accuracy, the spheres must be separated by at 
least 1 meter and be rotating at no more than 10 radians per second. Since the rotation causes a 
sinusoidal variation in the output, potential readings must be taken at least every 0. 1 seconds to get 
the desired resolution. Four booms (a 2.75m Weitzmann Quadrupole Stacer boom system) lie in a 
plane perpendicular to the spin axis. The fifth (a 1.65m Weitzmann Monopole Stacer boom) lies 
along the spin axis. A 5cm radius gold plated aluminum sphere is mounted at the end of each 
boom. 

3.3 Attitude Control and Determination System. The satellite is modelled as a spinning oblate 
platform which spins at 10 rpm. The orientation or attitude of each satellite will not be controlled, 
but the spin rate will be controlled to 10 ± 2rpm. An on-board cold gas propellant system will be 
used to spin-up the satellite initially as well as make any necessary spin rate adjustments. Two 2- 
axis magnetometers, 2 sun sensors, and 1 horizon crossing sensor are used to determine the 
attitude of each ELF satellite. By using different combinations of these 5 sensors, the attitude of 
the satellite can be determined at all times during the orbit Some of the data will be redundant but 
this redundancy can be used to enhance the accuracy of the readings. These readings will lie 
within the ± 1 degree error in attitude knowledge margin. 

3.4 Data Processing System. The data processing system is sized to store up to 24 hours worth of 
data. These data include the electrical field potential as well as the attitude readings. This system 
will also handle housekeeping functions on-board the satellite. Sampling 6, 12 bit words every 
0.1 seconds leads to a daily on board storage capability of only 7.8 M bytes. 

3.5 Communication System. Because the attitude of each ELF satellite will not be controlled and 
each satellite can maintain a different orientation, a virtually omni-directional antenna is needed on- 
board for communication. A stripline wraparound antenna meets this requirement Frequencies in 
the S-band will be used for receiving instructions and transmitting the collected data to a ground 
station. Data will be transmitted twice per day to one ground station. The ground station will use a 
4.3 m parabolic dish with tracking capabilities. The actual location of this ground station is yet to 
be determined. Each satellite will take a maximum of 82 seconds to transmit the stored data. It is 
expected that each satellite will pass within range of the ground station at least twice per day. 

3.6 Power System. The satellite’s power will come from solar cells wrapped around the exterior 
of the spacecraft. Since the satellites will not generally be placed in sun synchronous orbits, they 
will have to function in the dark as well. Therefore, the solar cells will be backed up with 
batteries. The minimum power generation will be 12.77 W which will be sufficient to cover the 
power requirements of all systems. 

3.7 ELF Satellite Structure and Configuration. The cylindrical primary structure is 45 cm in 
diameter and 35 cm high. It will be composed of 0.16 to 0.32 cm thick aluminum honeycomb. 
Most subsystem components will be mounted on this plate as shown by Figure 4. Individual 



component covers will provide radiation shielding as required. 

A 50 ELF SATELLITE MISSION 

Having designed the ELF satellite, the next step is to consider the ELF constellation deployment. 
Two different scenarios were considered by Redd and Olsen [1990]. The first considers using 
Orbital Sciences Corporation's Pegasus as a dedicated launch vehicle. The second option looks at 
piggyback opportunities on McDonnell Douglas's Delta II launch vehicle. The following 
discussion centers on the Pegasus vehicle to provide 8 launches with 6 satellites per launch. 
Pegasus will insert the payloads into 550km altitu de orbits for a minimum lifetime orbit of 7 years. 
When die proper orbit is reached, one of the 6 satellites will be ejected. The rest of the cluster will 
remain attached to the upper stage. This stage will then be manouvered to a slightly different 
altitude ( about 75km higher) or a slightly different inclination (about 2.5 degrees) where another 
ELF will be ejected. This sequence wifi be repeated until all ELFs are deployed. By slightly 
changing the altitude and/or inclination of each satellite, the satellites will be dispersed further by 
orbital perturbations caused by the earth's oblateness. 

Figure 5a shows the initial 8 orbit planes viewed from above the north pole. After six months each 
set of 6 satellites are dispersed from these orbits to produce the coverage shown in Figure 5b. In 
order to make the improved scientific measurements data from all satellites will be compared during 
a small time interval. Figure 5c shows the resulting satellite orbit segments for a 20 minute time 
interval. The ensuing observational coverage is the key parameter in deciding how to select the 
initial 8 Pegasus launches. 

CONCLUSION 

The concept of globally "imaging" the magnetospheric electric field at ionospheric altitudes (500 to 
1000km) has been considered. A constellation of 50 small ELF satellites simultaneously measuring 
the in-situ electric field has been outlined. This small ELF satellite has been designed to minimize 
mass and cost (Redd and Olsen, 1990). Each ELF satellite has a mass of 19 kg and costs less than 
$250,000 for the hardware. These ELF satellites would be lauched six at a time by Pegasus or four 
at a time as a secondary payload by the Delta n. Based upon the Pegasus launch senario, where 
one vehicle costs $7 million, a total of 9 launches would be needed. This would result in a mission 
cost of $90 million (assuming each ELF costs $0.5 million) or $171 million (assuming each ELF 
costs $2.0 million). Such a mission would fall in die NASA small to moderate mission category 
while the scientific data from such a mission would have a major impact on the ionospheric, 
thermospheric, and magnetospheric disiplines. 

Acknowledgments : This research was supported by NASA contracts NAGW-77 and NAGW- 
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Figure 3 Close-up exterior view of an ELF satellite 




Figure 4 Internal layout of the subsystem components 








Figure 5a View from the North Pole of the orbits after all 
8 Pegasus launches. Satellites are still in clusters of 6. 



Figure 5c Global coverage during a 20 minute time span 
6 months after deployment 
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ABSTRACT 


Dynamics Explorer (DE)-2 electric field and particle data have been used to constrain the 
inputs of a Time-Dependent Ionospheric Model (TDIM) for a simulation of the ionosphere on 22 
November 1981. The simulated densities have then been critically compared with the DE-2 electron 
density observations. This comparison uncovers a model-data disagreement in the morning sector 
trough; generally good agreement of the background density in the polar cap and evening sector 
trough; and a difficulty in modelling the observed polar F-layer patches. From this comparison, 
the consequences of structure in the electric field and precipitation inputs can be seen. This is 
further highlighted during a substorm period for which DE-1 auroral images were available. Using 
these images, a revised dynamic particle precipitation pattern was used in the ionospheric model; 
the resulting densities were different from the original simulation. With this revised dynamic 
precipitation model, improved density agreement is obtained in the auroral/polar regions where the 
plasma convection is not stagnant However, the dynamic study also reveals a difficulty of 
matching dynamic auroral patterns with static empirical convection patterns. In this case, the 
matching of the models produced intense auroral precipitation in a stagnation region, which, in 
turn, led to exceedingly large TDIM densities. 
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L INTRODUCTION 

Theoretical models of the terrestrial ionosphere and thermosphere have now attained global- 
scale [Sojka and Schunk, 1985; Roble et al., 1987; Fuller-Rowell et al., 1987]. Such models 
require global-scale inputs, especially for the magnetospheric parameters. At the present time, no 
global-scale theoretical model of the magnetosphere exists. Hence, statistical (empirical) models of 
the magnetospheric electric field [Volland, 1978; Heelis et al., 1982; Sojka et al., 1986; Heppner 
and Maynard, 1987] and auroral precipitation [Spiro etal., 1982; Hardy et al, 1985; Evans, 1987] 
are used as inputs to ionospheric and thermospheric models. Such studies are handicapped 
because statistical inputs cannot readily be adapted to simulate the spatial variations and temporal 
dynamics of geomagnetic storms and substorms. To minimize some of these difficulties, models 
have used in-situ satellite observations of auroral precipitation and electric fields, combined with 
Defense Meteorological Satellite Program (DMSP) images and ground-based incoherent scatter 
radar observations of electric fields, to define the auroral and convection boundaries [Sojka et al., 
1983; Rasmussen et al., 1986b], With the availability of global auroral images from the Dynamics 
Explorer 1 (DE-1) satellite, it is possible to directly define a magnetospheric auroral input on the 
scale needed by the ionospheric and thermospheric models. 

In a previous study, we investigated methods to provide a description of the energy flux 
and characteristic energy of the auroral precipitation based on DE-1 images [Sojka et al., 1989]. 
The resulting precipitation models were compared with in-situ DE-2 particle measurements. Rees 
et al. [1988] have carried out image-to-precipitation calculations and used similar in-situ cross 
checks on the auroral precipitation. In our previous study, a particularly dynamic auroral period 
was chosen because it provided excellent image sensitivity. Unfortunately, the magnetospheric 
convection was also dynamic, and because of the northward turning of the IMF during the selected 
period, a realistic convection model could not be defined. The study did, however, show the 
ionospheric dependence on the auroral dynamics. 

For this study we used a period where both the precipitation and convection patterns could 
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be compared with DE data. In addition, the IMF data was available and was found to be 
reasonably behaved. 

Section 2 describes the selected data period, while Section 3 describes the quiet day inputs 
and model-observation comparison. In Section 4, the major substorm disturbance in this period is 
analyzed. Results of this study are summarized in Section 5. 
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n. 22 NOVEMBER 1981 

This day was selected for a variety of reasons based upon both data availability and 

geomagnetic conditions. During this 24-hour period, the three-hourly Kp dropped from 3+ to 2 

(see top panel in Figure 1); the average three-hourly Kp was 2+- Several substorms occurred 
during this day as shown by the short duration peaks in the AE index in the bottom panel of Figure 
1. Between 0730 and 0935 UT at least two substorms occurred and the AE reached values of 750 
during this time. This was the strongest disturbance; for the rest of the period the AE index 
averaged around 200 with brief departures to nearly 500. These departures being associated with 
other smaller substorms. Because these substorms were only moderate intensity they would not 
lead to a significant perturbation of the auroral F-region. At 0150, 1540, and 2300 UT the DE-1 
auroral images show the substorms and confirm they are moderate compared to the 0730 and 0900 
UT substorms. Hence, the 0730 to the 0935 UT period will be considered in detail as a substorm 
period while the other substorms will be averaged over and considered as a quiet period. 

The Interplanetary Magnetic Field (IMF) data were also available and two components are 
shown in Figure 1 (middle panel). For almost the entire period the IMF had the same sector 
morphology, namely southward IMF with By positive. In Figure 1 (middle panel), the By 
component is about 4 nT (gamma) positive until about 2100 UT when it approaches zero and then 
it oscillates about this value for the remaining three hours. The Bi component is negative 
(southward) all day with the exception of a one-hour period from 1030 to 1 130 UT. During this 
brief period, the Bz component is weakly positive, less than one nT. There appears to be no 
correlation between the northward turning at ~1 130 and subsequent southward turning of the IMF 
at 1330 UT and the substorms, as indicated by the AE index. However, the decreasing Kp is 
probably correlated with the Bz component approaching zero from its initial -5 nT value. 

The IMF southward condition implies that a two-cell convection pattern is appropriate, 
while the long period of By positive implies that the convection pattern would have a well defined 
skewing. This corresponds to the Heppner and Maynard [1987] BC pattern. In our previous DE- 
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TDIM study [Sojka et al„ 1989], it was not possible to define the convection pattern. However, 
this study begins with the convection pattern established; this will be demonstrated in the next 
section. 

A further reason for selecting this period was the availability of DE-1 auroral images from 
the S AI instrument The shaded bars on the UT time axis of Figure 1 indicate the four periods for 
which auroral images were taken. Our attention will be focused upon the 0635 to 0935 UT period 
when substorms were present A third reason for focusing on this day was that it can be well 
represented by a "quiet" day with a three-hour substorm period superimposed. Finally, during this 
24-hour period data were available from a number of DE-2 orbits for almost all the instruments. 
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m. QUIET DAY STUDY 

The 22nd of November 1981 can be characterized as a quiet day with a superimposed pair 
of substorms occurring at 0730 and 0900 UT. In this section, the quiet day aspect of this 24-hour 
period will be analyzed. Data from several instruments on DE-2 will be used to help define the 
magnetospheric drivers for the ionosphere and, subsequently, to compare the model ionospheric 
densities with the in-situ observed densities. During this 24-hour period, data is available from 1 1 
partial DE-2 orbits out of a possible 15 orbits. Five of these northern hemisphere passes have been 
selected for detailed presentation because complete data sets for the Vector Electric Field Instrument 
(VEFI), Low Altitude Plasma Instrument (LAPI), and Langmuir Probe (LANG) are available. 
These five orbits represent the variability found during this 24-hour period, including one pass 
during the substorm . Figure 2 shows the selected passes, approximated as straight lines, in the 
polar magnetic coordinates with a Hardy et al. [1985] auroral oval shown as the shaded region. 
Each DE-2 passage is labeled with the approximate time at which the orbit crossed the dawn-side 
northern high latitude region. This time will be used as the reference for discussions about specific 
orbits. 

3.1. Electric Field 

Data from the VEFI instrument was used to verify the use of the Heppner and Maynard 
[1987] BC electric field model to simulate plasma convection on this day. The DE-2 ac and dc 
electric field instrumentation [Maynard et al., 1981] was designed for three-axis measurements 
using the symmetric double-probe floating potential technique with cylindrical antennas providing 
21.4 m baselines. As described by Maynard et al. [1982], the Z-axis (east-west geographic) 
antenna failed to deploy; hence, measurements were restricted to the two components in the orbit 
plane. The dc data for most purposes are presented in spacecraft coordinates in which the two 

components are resolved into a horizontal component Ex (positive in the direction of the spacecraft 
velocity vector) and a vertically upward component. A correction for contact potentials is 
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determined for each pass using an assumption that the 45° latitude points are at zero potential. For 

presentation and analysis, the Ex component is examined after v x B fields from the satellite's 
motion and co-rotation have been removed. 

From Figure 2 it can be argued that the missing electric field component on this day is of 
secondary importance because the orbits lie primarily along the dawn-dusk meridian. The 

dominant electric field component in this region is along the VEFI Ex direction. This electric field 
component is available and, hence, its path integral gives a reasonable approximation to the cross 
polar-cap potential. For this day there are no DE-2 retarding potential analyzer (RPA) data to infer 
the missing electric field component. 

Five passes of VEFI data were compared with the Heppner and Maynard [1987] BC 
electric field model (HM-BC). Figure 3 shows the comparison between the VEFI and HM-BC 
electric potential for the five orbits shown in Figure 2. In this figure, each panel corresponds to a 
single pass and is labeled with the orbit time. The HM-BC pattern has the strong electric field in 
the dawn sector of the polar cap; this asymmetry is associated with the positive IMF By . In the 
potential plots shown in Figure 3, the strong electric fields are associated with regions of large 
gradients in the potential. The VEFI data, especially the top two panels, show even larger 
asymmetries in the electric field. This is consistent with the IMF data (Figure 1), which shows that 
the By component is strongest for the first 8 hours of this day. Toward the end of this period, later 
than 2100 UT, the IMF Bi and By components become very weak and the convection pattern is 
expected to depart from the BC pattern. This change is also manifested in the decrease in Kp 
(Figure 1). The Heppner and Maynard potential patterns have been parameterized and made Kp 
dependent by Rich and Maynard [ 1989], but we have not included this variation in our simulation. 
In this instance we believe that the departure from the BC convection pattern as Bz approaches zero 
is more significant than the 30% reduction in Kp, and because we don't know how the convection 
pattern changed, we adopt a constant BC pattern for the quiet period. The comparison of the 
electric potential magnitude is best at the latter times when the IQ) is lower. At the earlier times, the 
electric potential variations observed by VEFI are larger than the BC model. By simply changing 
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the Kp better agreement is not obtained, because the model smooths out these regions of strong 
electric fields. 

3.2. Auroral Precipitation 

The Hardy et al. [1985] auroral precipitation model for Kp=3 was compared to the DE-2 
plasma observations. The LAPI instrument on DE-2 consisted of 15 divergent parabolic-plate 
electrostatic analyzers, each of which provided differential spectral measurements of electrons and 
positive ions over the energy range from ~ 5 eV to ~ 30 keV with a 1 sec temporal resolution. The 
analyzer ficlds-of-view were selected to provide optimal sampling within and outside of the particle 
source and loss cones at DE-2 altitudes (500—1000 km). For the observations on 22 November 
1981, 31-point energy spectral measurements were made for electrons at pitch angles 0°, 7.5°, 15°, 
30°, 45°, 60°, 105°, 135°, 165°, and 172.5°. Further details and calibration information on LAPI 
may be found in Winningham et al. [1981]. 

For each LAPI spectrum, the precipitating electron energy flux from 50 to 20,000 eV was 
integrated to give an auroral energy flux. Figure 4 shows these energy fluxes for the five selected 
orbits. Also shown in this figure is the Hardy et al. [1985], Kp=3 energy flux at the satellite 
locations. The auroral precipitation is both more structured and more confined than the statistical 
oval. For the first half of the day, the LAPI observations are comparable to the model predictions 
in the oval regions. However, toward the end of the day as the activity decreases (see Figure 1), 
the LAPI data show lower fluxes of precipitating electrons than the model. In the first three passes 
shown in Figure 4, the model fluxes are in conflict with the observations in the polar region, i.e., 

no precipitation is observed while the model suggests about 0.2 erg/cm 2 s should be present. This 
descrepancy is not severe when it is realized that these three passes are not truly in the polar cap, 
but just touch the dayside cusp region (see Figure 2). In this region, the difficulty of using a 
statistical model is highlighted by the model's dayside oval width being significantly wider in 
latitude than was actually present The remaining two orbits do cross the polar cap and the fluxes 
do decrease as expected (see Figure 4, bottom two panels). Hence, the polar cap problem in the 
first three passes in Figure 4 indicates that the actual oval was somewhat more expanded on the 
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dayside than the Kp=3 oval we have adopted. This is consistent with the magnetic indices, which 

were initially at a Kp=3+ level. For the purposes of this study, the Kp=3 oval appears to be 
reasonable and, hence, we adopted the Hardy Kp=3 model to describe the auroral precipitation for 
the quiet day and for consistency with the adoption of a fixed convection pattern. 

3.3. Ionospheric Model 

Our ionospheric model was initially developed as a mid-latitude, multi-ion (NO+, 02 + , 

N2+, and 0+) model by Schunk and Walker [1973]. The time-dependent ion continuity and 
momentum equations were solved as a function of altitude for a co-rotating plasma flux tube 
including diurnal variations and all relevant E and F region processes. This model was extended to 
include high latitude effects due to convection electric fields and particle precipitation by Schunk et 
al. [1975, 1976]. A simplified ion energy equation was also added which was based on the 
assumption that local heating and cooling processes dominate (valid below 500 km). The 
trajectories for plasma motion were determined as the plasma moved in response to convection 

electric fields. A further extension of the model to include the minor ions N+ and He+, an updated 
photochemical scheme, and the mass spectrometer/incoherent scatter (MSIS) atmospheric model is 
described in Schunk and Raitt [1980]. 

The addition of "constant" plasma convection and empirical models for particle precipitation 
is described in Sojka et al. [1981a, b ]. More recently, the ionospheric model has been extended 
by Schunk and Sojka [1982] to include ion thermal conduction and diffusion-thermal heat flow so 
that the ion temperature is now rigorously calculated at all altitudes between 120 and 1000 km. 
The adopted ion energy equation and conductivities are those given by Conrad and Schunk 
[1979]. Also, time-dependent plasma convection and particle precipitation inputs have been used 
with the basic high-latitude model so that magnetospheric storm and substorm conditions could be 
studied [Sojka and Schunk , 1983, 1984]. 

For the quiet day aspect of this study, the above ionospheric model was run for a 24-hour 
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period during which the magnetospheric and solar inputs were fixed. This produced a "diumally 
reproducible" electron density data set. The solar conditions were represented by a solar flux 
F10.7 value of 180. Based upon the indices in Figure 1 and the DE data shown in the previous 
two sections, the magnetospheric inputs (convection and precipitation) were, respectively, the 
Heppner and Maynard [1987] BC model and the Hardy et al. [1985] Kp=3 model. Since both of 
these models are extensively discussed in the respective references, no further discussion is given 
here. A neutral atmosphere is also needed; the MSIS model was used for this input [Hedin, 1987]. 
For the MSIS model, the above FI 0.7 value and Ap=l 1 were used as input parameters. A simple 
meridional 200 m/s neutral wind model blowing from 1300 MLT to 0100 MLT was adopted for 
the dark polar regions. The validity of this wind is difficult to quantify since no IMF dependent 
wind model is presently available. Rees and Fuller-Rowell [1988] used their coupled 
thermosphere and ionosphere to study the IMF By dependence of the high latitude. At latitudes 
poleward of ~75\ a strong asymmetry is found in the F-region neutral wind as a function of IMF 
By, but at latitudes equatorward of this a meridional wind peaking at ~0100 MLT is found. Since 
the wind only affects the F-layer at the latitudes equatorward of 75*, this simple wind is consistent 
with that from the more physical coupled ionosphere-thermosphere model. The TDIM produced a 
data base of electron densities for magnetic latitudes poleward of 50* over the 100 to 800 km 
altitude range. 

The quiet-day ionospheric simulation for 22 November 1981 is shown in Plate 1 as a UT 
series of snapshots of NmF2 (top panels), hmFl ( 2nd from top panels), molecular ion density at 

300 km (third row), and 0+ at 800 km (bottom panels). Six UT snapshots uniformly distributed 
throughout the day are shown. N m F2 indicates how significant the UT control of the polar 

ionosphere is: at 0100 UT the polar densities are below 4 x 10 5 cm-3, while at 1700 UT the polar 

densities everywhere approach 1 X 10^ cm-3. Apart from the absolute density differences, the 
polar cap density is structured. A "tongue of ionization" is evident at most times and is seen as the 
channel of high densities moving through the cusp region into the polar cap. This tongue extends 
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first in a westward direction and then over the pole in a predominantly anti-sunward direction. The 
dawn polar cap has significantly higher densities than does the dusk polar cap. All these density 
structures are associated with the convection pattern, the IMF By positive Heppner and Maynard 
pattern. The density variation with UT is controlled by the location of the transition from daylight 
produced ionization to darkness (solar terminator) in relation to the cusp region where plasma 
enters the polar cap. As this terminator region approaches the cusp, the density in the polar cap 
increases (compare densities at 0500 and 1700 UT). This effect has been observed by ground 
-based incoherent scatter radars [de la Beaujardiere et ai, 1985]. 

The second panel from the top shows the height of the F-region peak density ( hmF2 ). In 
sunlight, the peak lies just below 300 km, whereas inside the dark polar cap it ranges from as low 
as 220 km up to 400 km. The highest values are located in the polar cap dawn sector just poleward 
of the cusp. In this region, the electric fields are the largest (largest flow speeds) and the flow is 
predominantly towards the magnetic pole. This results in a strong upward drift, i.e., the electric 
field is perpendicular to the tilted magnetic field which gives rise to an upward E x B drift. 
However, on the other side of the pole, this flow leads to a downward drift which appears as the 

region of lowest hmF2 values. 

In the vicinity of the peak (300 km), the molecular ions do not show significant polar cap 
densities or structure (third row). The molecular ions instead, are enhanced in the auroral oval and 
the sunlit region. This is consistent with the fast recombination time for the molecular ions which 
inhibits transport effects. Only a very small transport enhancement in the molecular ions can be 
seen in the region of the high hmF2 where the plasma moves into the polar cap. 

At 800 km, the ion density variation (bottom panel, Plate 1) is similar to the NmF2 vari- 
ation. For a diumally reproducible study such as this quiet day, this is not surprising since the 
topside time constant, on the order of an hour, is 'fast' compared to the times at which the model 
inputs are changing. This would not be the case during a substorm [Sojka and Schunk , 1983, 
1984]. 
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3.4. Electron Density Comparison 

The computed electron densities were compared with those measured by the Langmuir 
Probe (LANG) on DE-2. The Langmuir Probe instrument was designed to perform in-situ 
measurements of the electron temperature and the electron and ion densities in the ionosphere. 
Two independent sensors were connected to individual adaptive sweep voltage circuits and, hence, 
the LP can continuously track the changing electron temperature and spacecraft potential while 
autoranging electrometers adjust their gain in response to the changing plasma density. The control 
signals used to achieve this automatic tracking provide a continuous monitor of the ionospheric 
parameters. Additionally, internal data storage circuits permit high resolution, high data rate 
sampling of selected volt-ampere curves for transmission to the ground to verify or correct the 
inflight processed data. The availability of analog telemetry channels in the DE spacecraft permits 
the transmission of raw electrometer data from either of the two sensors. Further details can be 
found in Krehbiel et al. [198 1]. 

Figure 5 shows the electron density along the five orbits for the LANG and the TDIM. In 

each panel the logarithm to the base 10 of the density (cm-3) is plotted against UT. Overall, the 
model densities track the observed densities fairly well, but there are a few exceptions to this. In 
the dawn sector at the beginning of each pass, the model is consistently lower than the 
observations (see arrows labeled A). The LANG density structure in the center of the polar cap is 
also missing from the model (see arrows labeled B) and the location of the structure in the dusk 
sector trough (arrows labeled C) does not always coincide with the observed location. Ignoring 
for the time being the dawn trough discrepancy (A), the density comparison for these five orbits 
ranges from good to almost excellent (orbit 1400 UT). This is also the case for the other orbits 
(not shown) for which LANG data were available where the differences fall into the same 
categories found in Figure 5. These differences are discussed and related to the geophysical 
conditions in the following section. 
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3.5 Quiet-Day Discussion 

Both the LANG observations and the model show that the density does not peak in the 
auroral region, but rather the peak density is found in the polar cap (see Figure 5). This is because 
the densities are all in the topside ionosphere at the satellite altitude where the role of plasma 
transport is more important than that of the auroral precipitation [Sojka and Schunk , 1984]. The 
DE-2 altitude ranges from 800 km on the dawnside to 400 km on the dusk end of each passage 
(see Figure 2), and hence, the measurements were made above the F-region peak. Agreement is 
also present in the dynamic range of the observed densities with, of course, the exception of the 
dawn trough region. 

The disagreement in observed and modelled densities found in the morning sector trough 
region was unexpected, since this region is believed to be relatively independent of the magneto- 
spheric inputs. Figure 6 shows the relevant ionospheric morphology in this region. A circle 
represents the region where the model densities are low relative to the LANG densities. This circle 
is located equatorward of the nocturnal auroral region and is in darkness; the straight line 
represents the solar terminator. Passing through this region is the 0735 UT DE-2 satellite track. In 
the TDIM, the plasma present in the large circle has co-rotated to this region from around midnight 
During this co-rotation the plasma flux tubes have remained equatorward of the auroral oval (see 
dashed line in Figure 6). These co-rotational and auroral conditions are consistent with the VEFI 
observations (Figure 3) in that the electric potential is not changing with latitude (E * 0) and the 
LAPI (Figure 4) data which indicates the region to be equatorward of the auroral precipitation. The 
TDIM densities in these flux tubes are maintained by the upward induced drift caused by the 
generally equatorward neural winds. Although this mechanism is adequate to maintain the 
ionosphere around midnight, it is unable to do so just prior to sunrise [Sojka et al ., 1981b]. This 
is the first attempt to check the topside densities in the trough at these pre-sunrise times. The 
model densities are clearly too low in this region for the first four orbits shown in Figure 5, 
although the final orbit does agree better with the model. An extensive check has been made to 
confirm that the observed densities are accurate, and LANG data from other days confirms this 
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trend. This result has led us to begin a follow-up study to determine what is missing in the TDIM 
nighttime trough processes. The obvious candidate is a downward plasma flux from the 
plasmasphere, which was set to zero in the simulation owing to a lack of measurements of this 
parameter. 

Polar cap density structures, as highlighted by the B arrows in Figure 5, are not expected to 
be present in the TDIM. The model uses only statistically smooth inputs for the magnetospheric 
convection and precipitation and, hence, no source for multiple structures in the polar cap is 
present. However, this picture is somewhat over simplified; certain aspects of the structure can be 
related to the presently used inputs. For the first two orbits in Figure 5, the first B arrow points to 

a region in the polar cap where the density increases from about 3 x 10* to 3 x 105 cm-3 over a 
very small distance. Figure 3 VEFI data show that this is the region of largest electric fields in the 
polar cap. These measured electric fields are stronger than those obtained from the HM-BC input 
to the TDIM simulation. These electric fields are responsible for convecting plasma through the 
cusp region into the polar cap. The polar cap density features then depend upon the plasma and 

convection in the vicinity of the cusp. Such a cusp control can be seen in NmF2 as a function of 
UT in Plate 1. These convecting density structures have been reported [Weber et al., 1984, 1986] 
and are referred to as F-layer patches. Morphologically, they are present during IMF southward 
conditions, while their sizes and locations depend on the variability in the cusp convection in ways 
not fully understood. Earlier TDIM ionospheric studies have shown that significant polar cap 
density structuring is possible even with smooth inputs [Sojka and Schunk, 1987; Rasmussen et 
al., 1986a]. In this study, the HM-BC model does not have a strong enough electric field to 
reproduce the structures observed along the first two orbits shown in Figure 5. Density structure 
B' in the 2210 UT pass (Figure 5, bottom panel) is a larger feature. Towards the end of the quiet 
period, the IMF Bz component approached zero and, hence, the actual convection pattern must 
have changed. This density feature may well be the equivalent of the TDIM density enhancement 
from 2220 to 2223 UT in the same orbit panel of Figure 5. The B' density structure is certainly 
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not a precipitation feature (Figure 4) since no appreciable precipitation is present 

The C arrows in Figure 5 highlight the third region where differences are present between 
the modelled and observed densities. In this case, the differences are not systematic and are not 
present in all cases. This region is the afternoon sector trough. Here, the depleted trough densities 
are associated with a combination of the convection-co-rotation balance and the presence of the 
auroral equatorward boundary [Sojka et al ., 1981a]. For the orbit at 1400 UT, an almost perfect 
agreement is found in the densities (Figure 5) and the convection (Figure 3), while for the orbit at 
2210 UT, the observed trough is significantly more poleward and a similar difference is seen in the 
convection at that time. Hence, the differences highlighted by the C arrows are probably 
associated with the electric field variability not present in the HM-BC convection model. At 
present, there are no methods available to reconstruct the variability in the electric field data. 

What about the "substorms"? In Figure 5 the 0735 UT pass occurred during the substorm 
period. However, it does not look especially different from the other orbits. Neither the 
convection (Figure 3) nor the precipitation (Figure 4) appear significantly more disturbed. 
Nevertheless, the DE-1 spin-scan auroral imager (SAI) data do show significant differences in the 
precipitation during this time. The next section looks at this substorm period in depth. 
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IV. THE 0635-0935 UT SUBSTORM PERIOD 


4.1. Auroral Imagery 

The high-altitude satellite DE-1 obtained continuous sequences of auroral images in time 
intervals in excess of four hours for some orbits during the fall of 1981. The satellite is equipped 
with a spin-scan auroral imager (SAI) that comprises three scanning photometers with off-axis 
parabolic reflectors as the primary optical elements. Two of these imaging photometers provide 
global auroral images at visible wavelengths, each utilizing any of 12 filters selected by ground 
command. The third photometer is capable of auroral imaging at vacuum ultraviolet (VUV) 

wavelengths with a similar selection of one of its 12 filters. In typical operation, a 30° x 120° 
image is obtained every 12 minutes for each photometer. The images are centered in the Earth 
nadir direction, with the 30° angular width of each image aligned transverse to the plane of satellite 
rotation and the coplanar satellite orbital plane. A detailed description of the imaging 
instrumentation is given by Frank et al. [1981]. 

During the substorm period, images were obtained every 12 minutes from 0635 to 0912 
UT. A total of 14 images show the auroral dynamics during this period. Images were taken 
simultaneously at 557 and 630 nm in the visible spectrum and in the vacuum ultraviolet (VUV) 
over the 123—160 nm bandwidth. On this day the 630 nm intensities were very low and because 
of statistical considerations were not suitable for use in the energy flux-characteristic energy 
analysis described by Sojka et al. [1989]. Plate 2 shows the VUV images for this substorm period 
and the color key gives the observed SAI photometer compressed pixel counts. From these 
observations, the auroral dynamics is evident. In the evening and night sectors, discrete auroral 
forms are present over a large region and they display a significant temporal variation, i.e., 
changes in these structures are seen from one image to the next Substorm onsets are observed in 
the 0723 and then the 0836 UT images. None of these spatial or temporal features are present in 
the statistical auroral model used in the quiet-day study. 

To use the images in Plate 2 as an input to the TDIM, they must be converted to 
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precipitating energy fluxes. The lower panels in each pair of images in Plate 2 show these energy 
fluxes. Sojka et al. [1989] discussed at length the procedures used to convert the SAI images to 
energy flux. In this study, two levels of conversion were considered; the first is the relative energy 
flux which conserves the auroral spatial and temporal variability, and the second is the absolute 
energy flux calibration. As has been noted by Sojka et al. [1989] and others [Rees et al. 1988], 
the second part is extremely difficult to do from first principles. The procedure used here is to 
obtain the relative auroral dynamics (spatial and temporal) from the images and then to scale the 
energy fluxes to achieve consistency with the LAPI data. This follows the general procedures 
described by Sojka et al. [1989]. The energy fluxes were computed from the VUV images by the 
algorithm described by Rees et al. [1988]. These energy fluxes were within a factor of two of the 
LAPI energy fluxes, which statistically was the same as found by Sojka et al. [1989]. Hence, the 
same scaling factor was used to scale all the computed image energy fluxes to the LAPI level. The 
resulting energy fluxes are shown in the lower panel of each pair in Plate 2 with a color key in 

ergs/cm 2 s. The 14 images span the entire image period and during this time the auroral region 
undergoes significant changes. The AE index peaks at 0800 and again at 0930 UT, which roughly 
correspond to the times of maximum expansion in the substorms. 

4.2. Substorm Density Comparison 

The TDIM was rerun using energy fluxes obtained from DE-1 images to simulate the 
ionospheric response to the auroral dynamics during the substorm period. From Figure 1 it can be 
argued that prior to about 0635 UT the auroral conditions were quiet, i.e., the AE index was 
relatively small and varied smoothly. On this basis, the previous quiet-day TDIM densities were 
used as the initial conditions at 0635 UT for the substorm simulation. During the substorm period, 
insufficient data were available to define a dynamic convection model and, hence, the HM-BC 
model was used again. Therefore, only the auroral inputs were changed for the substorm 
simulation. The Hardy empirical auroral model was replaced with a time varying sequence of 
image energy fluxes. These images occur at 12 minute intervals. In using the image energy 
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fluxes, each was assigned an image UT (the middle of its 12 minute period) and then energy fluxes 
at a given location were obtained by linearly interpolating between two images that occurred before 
and after the required TDIM UT. 

In creating the dynamic auroral model from the SAI images, it was difficult to infer the 
characteristic energy of the precipitation on this day. This is because the characteristic energy is 
deduced from the ratio of emissions at different wavelengths. The 630 nm emission is crucial for 
this ratio calculation. As already pointed out, the 630 nm had low count rates, hence, these data 
could not constrain this parameter. The computed ionospheric densities do, however, depend 
upon the characteristic energy. In the Appendix, a parametric study of this dependence is given. 
Based on the limited data available, we adopted a characteristic energy model; the SAI images, the 
LAPI data, and the Hardy et al. [1985] characteristic energies were used to define this model. The 
model defines the characteristic energy in the polar cap to be 1 KeV and equatorward of the oval to 
be 5 KeV. Inside the oval in the morning sector it is also 5 KeV, while in the afternoon-evening 
oval it is 3 KeV. In this model the oval is defined as regions where the energy flux exceeded 0.5 

erg/cm 2 s. 

Figure 7 compares the substorm TDIM (Adopted Eo line), the quiet day TDIM (Original 
Simulation line), and the LANG observations for the 0735 UT DE-2 pass. The model/measure- 
ment discrepancy in the dawn sector trough is, as expected, unchanged. In the polar cap, 
especially the dusk sector, the new model densities are in significantly better agreement with the 
LANG observations. However, the earlier similarity between model and data in the dusk sector 
trough is changed and now there is a marked difference in densities. The substorm TDIM 

(Adopted Eo) densities are almost an order of magnitude larger than the LANG trough densities, 
whereas the quiet day TDIM densities were only a factor of two larger than the LANG densities. 
This feature is significant and will be discussed in the following section. A second northern 
hemisphere DE-2 LANG pass occurred later in the substorm period. Figure 8 shows the LANG 
density (crosses), as well as the quiet day (solid line, Original Simulation) and substorm (solid 

line. Adopted Eo) simulation densities. This pass occurred between 0915 and 0935 UT, which 
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was toward the end of the image sequence and at the time of the second peak in the AE index (see 
Figure 1). Because several hours have elapsed since the beginning of the substorm, the simulation 
is independent of the initial quiet conditions. This is in contrast to the pass at 0735 UT which 
occurred near the beginning of the image sequence. In Figure 8, the substorm simulation shows 
significantly better agreement with the observed densities. It lacks the sharp polar cap density 
structures, but this is to be expected due to the lack of convection structure in the model. As with 
Figure 7, the biggest disagreement occurs in the dusk sector auroral oval/trough region. In this 
case, the trough minimum is matched reasonably well; however, the model predicts a huge density 
maximum in the oval at 0928 UT. The following section describes in depth how the discrepancy 
in Figure 7 arose. The argument to be presented also applies to the density difference shown in 
Figure 8. 

4.3. Substorm Discussion 

The density comparison in the polar regions indicates that the dynamic auroral TDIM 
simulation yields a better agreement with the data (see Figure 7). However, of concern is the 
larger discrepancy produced in the dusk sector trough. In order to understand this. Figure 9 is 
used to graphically display the plasma mechanisms leading to this greater difference. The circle 
indicates the region where the substorm TDIM produced the larger densities. This region lies in a 
region of auroral emissions, as indicated by the shaded region in Figure 9. In comparison, the 
original oval does not coincide with this circle (compare Figure 9 with Figures 6 or 2). Hence, the 
images show that during the substorm period, auroral precipitation extends considerably more 
equatorward in this local time sector than predicted by the statistical oval. In fact, this 
interpretation is consistent with the LAPI data which show strong precipitation equatorward of the 
statistical oval (see Figure 4, the 0735 UT pass at about 0753 UT). 

Also shown in Figure 9 are the trajectories of two plasma flux tubes that form this 
enhanced density region. On each trajectory the small circles mark hourly intervals. One of the 
trajectories (labeled a ) is for plasma convecting around from the dayside and stagnating in the 
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circle. This trajectory is responsible for the high TDIM densities because the plasma stagnates in a 
region of high auroral energy fluxes. However, such a scenario is extremely unlikely; stagnation 
occurs only in regions equatorward of the auroral oval where the magnetospheric electric field 
tends to cancel the co-rotational electric field. Hence, the inability to dynamically modify the 
convection pattern during the substorm has resulted in the "spurious" high TDIM densities. 
Furthermore, as seen in Figure 7, there is almost no dusk sector trough in this substorm 
simulation. From this it can be deduced that the magnetospheric electric field must have penetrated 
to lower latitudes. The single DE-2 substorm pass at 0735 UT in Figure 3 indicates that this is the 
case. The second plasma flux tube trajectory (labeled b ) shown in Figure 9 represents the typical 
poleward edge of the dusk sector trough. Here, tire plasma is convecting sunward inside the oval 
to the region of interest High TDIM densities are associated with this flux tube because it is in a 
region of high precipitation for several hours. The DE-2 orbit track at 0735 UT is superimposed 
on Figure 9 as the line with small crosses. 

V. SUMMARY 

In this study we have used an extensive DE data base to both constrain inputs to the TDIM 
and then check the simulated densities. This has been done for both a quiet period and a substorm 
period. The quiet period had a duration of about 24 hours, while the substorm duration was 3 
hours and occurred 6.5 hours into the quiet period. This is the most extensive high latitude TDIM- 
data comparison carried out to date. Aspects of the study can almost be regarded as looking at 
ionospheric "weather" rather than the usual "climatology." Such a capability can only arise when 
the inputs constraining the model are available with adequate resolution and coverage. Under ideal 
conditions the combined DE-1 and DE-2 data bases can almost achieve this resolution, as the study 
has shown. 

The study can be summarized by the following key findings: 

(1) The quiet-day study produced very good agreement between modelled and 
observed electron densities in the topside ionosphere with two important 
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exceptions. 

(a) Across the polar region the DE-2 LANG densities showed fine structure (F- 
layer patches) in addition to the overall regional density morphology. Such 
fine structure is associated with the "weather" variability in the 
magnetospheric inputs. These cannot be completely described by the 
present day observations. The full extent of this structure is, therefore, 
missing from the simulated densities. 

(b) A surprising discrepancy arose in the pre-sunrise mid-latitude trough. The 
TDIM densities were an order of magnitude lower than those observed by 
DE-2. This discrepancy was not restricted to this day; the measured 
densities were similar in this region over many days. This region is also 
independent of the magnetospheric inputs and so does not rely on good 
resolution in the DE inputs. A follow-up study has been initiated to elucidate 
the source of this marked discrepancy. 

(2) The substorm study showed remarkably good agreement with the observed 

densities with again two exceptions. 

(a) The observed fine structuring of the polar cap densities was even more 
marked during the substorm period. During a substorm, the highest degree 
of "weather" is present. No method of achieving both the temporal and 
spatial coverage is available to describe the source of this structure. 

(b) Although the auroral image sequence produces a very good description of 
the global substorm, the associated convection electric field is not available. 
This can lead to drastic errors in the model prediction. In this case, the dusk 
sector stagnation occurred inside the auroral oval. The ensuing TDIM 
density enhancement was over an order of magnitude greater than observed. 
Such a discrepancy is purely an artifact of not having the necessary electric 
field measurements to couple the auroral and convection patterns self- 
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consistently. 

This study shows again the importance of using extensive data bases to constrain the global 

and temporal behavior of the magnetospheric inputs to the ionospheric (and thermospheric) 

system. Even under quiet conditions, as shown by this study, simplifying model assumptions 

need to be verified. Under dynamic substorm conditions, even simple mismatching between inputs 

leads to major errors in calculated ionospheric densities. More work is needed in deducing an 

empirical means of making convection and precipitation patterns self-consistent during temporal 

variations. This would be even more useful as global auroral imagery becomes a standard 

replacement for statistical auroral inputs. The major "weather" aspects will require not only high 

resolution imagery of the precipitation, but also information concerning the structure in the 

convection pattern. It is the structure in the convection pattern that produces the density fine 

structure, as shown simplistically in the TDIM polar cap structure of Plate 1. 
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APPENDIX 

The depth into the thermosphere that an auroral electron can penetrate and the altitude 
distribution of ionization depend upon several factors. The most important factor is the energy of 
the electron. When a flux of electrons travels through the ionosphere, the energy is characterized 
not by a single value, but by the initial energy spectrum. A lengthly calculation is needed to deduce 
the associated altitude distribution of ionization. To simulate this process, the TDIM selects 
between ionization profiles for different characteristic energies, different seasons, and different 
solar cycle conditions. These basic ionization profiles were computed using the Strickland auroral 
deposition code (R. Daniel, private communication, 1988). 

Figure 10 shows the sensitivity of the TDIM density calculation to the characteristic energy. 
The substorm simulation described in Section 4 was repeated three times for three fixed charac- 
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teristic energies. In each case, the energy fluxes varied as shown in Plate 2, i.e., according to the 
image variability. Three results for typical auroral energies are shown in Figure 10, namely 1, 3, 
and 5 KeV. Since the densities shown in Figure 10 are along the DE-2 track used in Figure 7, they 
are all in the topside ionosphere. In this region, the lowest energies contribute the most ionization, 
although the integrated ionization is largest for the most energetic electrons. The curves shown in 
Figure 10 indicate that a factor of 2 to 4 density difference can be anticipated for the expected 
uncertainty in the auroral electron energy. 

For the Adopted Eo data shown in Figures 7 and 8, the following simple auroral charac- 
teristic energy distribution was adopted: 

1) Inside the polar cap, poleward of the auroral oval magnetic latitude, Eo = 1 KeV. 

2) Inside the afternoon-evening sector oval, Eo = 3 KeV. 

3) Inside the morning sector oval, Eo = 5 KeV. 

4) Equatorward of oval, Eo = 5 KeV. 

In this case the oval is defined as a region where the auroral energy flux exceeds 0.5 erg/cm2s. 
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ILLUSTRATIONS 


Figure 1. Geomagnetic and IMF indices for 22 November 1981, and a chart indicating when DE- 
1 auroral images were available (gray shading on the UT axis). The 3-hourly Kp index is plotted 
in the top panel while the auroral electroject index, AE, is plotted in the bottom panel. Two 
components of the IMF, By and Bz, are plotted in the middle panel. 

Figure 2. DE-2 orbit tracks in the northern hemisphere in dipole latitude-MLT coordinates in a 
polar format These five orbits are used to highlight the electric field and particle precipitation data. 
In the dawn sector, the DE-2 altitude is about 800 km which decreases to below 400 km in the 
dusk sector. A Hardy et al. [1985] Kp=3 oval is shown for reference, as the shaded region. The 

shading corresponds to energy fluxes larger than 0.25 ergs/cm 2 s. 

Figure 3. Electric potential along the DE-2 satellite track for five northern hemisphere passes on 
22 November 1981. The electric potential computed from the VEFI electric field data (dashed line) 
is compared with that from the Heppner and Maynard [1987] BC Kp=3 model (solid line) is 
shown as the solid line. The five satellite tracks are shown in Figure 2. 

Figure 4. Electron precipitation flux along the DE-2 satellite track for five northern hemisphere 
passes on 22 November 1981. The precipitation flux computed from the LAPI electron instrument 
(dashed line) is compared with that from the Hardy et al. [1985] Kp=3 model (solid line). The 

energy flux is in ergs/cm 2 s. The five satellite tracks are identical to those used in Figure 3 and 
shown in Figure 2. 

Figure 5. Electron density along the DE-2 satellite track for five northern hemisphere passes on 
/ 22 November 1981. The electron density from the LANG instrument (" + " symbols) is compared 
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with that from the TDIM (solid line). The electron density is in cm- 3 . The five satellite tracks arc 
identical to those in Figures 2-4. 

Figure 6. Schematic view of the plasma/DE-2 geometry in the low density pre-dawn trough region 
plotted in a dipole latitude-MLT polar format. This region is shown as a circle on the line repre- 
senting a DE-2 satellite track, and the dashed line represents a co-rotational plasma trajectory. The 
shaded region shows the Hardy et al. [1985] auroral oval, while the terminator is shown as the 
solid line across the dayside. The format is the same as for a single panel of Figure 5. 

Figure 7. Electron density comparison during the 0735 UT substorm on 22 November 1981. 
The DE-2 LANG data are compared with the TDIM simulations. Two TDIM simulations are 

shown, the original (from Figure 5) and the storm case labeled Adopted Eo. 

Figure 8. Second substorm electron density comparison on 22 November 1981. The DE-2 
LANG data are compared with the TDIM simulations. The figure layout is identical to that of 
Figure 7. 

Figure 9. Schematic view of the plasma/DE-2 geometry in the "stagnation" dusk trough region 
plotted in a dipole latitude-MLT polar format This region is shown as the intersection of the solid, 
crossed, line representing a DE-2 satellite track and the solid, circled line (labeled a ) representing a 
plasma trajectory. The shaded region shows the auroral precipitation inferred from an appropriate 

DE-1 SAI auroral image (> 0.25 ergs/cm 2 s). 

Figure 10. Model electron density variation along the DE-2 orbit track during the 0735 UT storm 
period. Each curve represents the TDIM simulation for different characteristic energies. In each 
case the SAI images were used to determine the energy flux for the auroral precipitation. The 
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densities are in electrons/cm3. This is the same orbit shown in Figure 7. 

Plate 1. Color coded snapshots of the quiet-day ionospheric simulation at six equally spaced 
times. Top row corresponds to NmFi, 2nd row to hmF2, 3rd row to the molecular ion density at 
300 km, and the bottom row to the electron density at 800 km. Data are color coded in magnetic 
latitude-MLT polar format from 50 degrees to the pole. The polar dial format is identical to that 
used in Figures 2, 6, and 9. Color keys are given on the right side of each row. 

Plate 2. Fourteen DE-1 VUV auroral images and the inferred auroral energy fluxes during the 
substorm period on 22 November 1981. The top panel in each pair is the VUV compressed count 
rate image, while the lower panel is the energy flux in the same polar plot used in Plate 1. Color 
keys are given on the right side of each row. 
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